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内 容 简介 


本 书 从 数学 建 模 的 角度 介绍 MATLAB 的 应 用 。 本 书 的 作者 都 具有 实际 的 数学 建 模 参 赛 经 历 
和 竞赛 指导 经 验 。 书 中 内 容 完全 是 根据 数学 建 模 竞 赛 的 需要 而 编排 的 ,涵盖 了 绝 大 部 分 数学 建 模 
问题 的 MATLAB 求 解 方法 。 全 书 内 容 分 上 下 两 篇 。 上 篇 介绍 数学 建 模 中 常规 方法 的 MATLAB 
实现 ,包括 MATLAB 交互 ,数据 建 模 、 程 序 绘图 \ 灰 色 预 测 、 规 划 模型 等 方法 ;还 介绍 了 各 种 高 级 方 
法 的 MATLAB 实现 ,包括 遗传 算法 、 粒 子 群 算法 、 模 拟 退火 算法 人工 神 经 网 络 ,小 波 分 析 、\ 动 态 仿 
真 , 数 值 模拟 等 。 下 简 以 真实 的 数学 建 模 赛 是 为 案例 ,介绍 了 如 何 用 MATLAB 求 解 实际 的 数学 建 
模 问题 ,给 出 了 详细 的 建 模 过 程 和 程序 。 书 中 的 附件 部 分 介绍 了 作者 在 建 模 竞赛 中 屡 获 大 奖 的 经 


验 。 相 信 这 些 经 验 对 准备 参加 数学 建 模 竞赛 的 读者 会 有 所 帮助 。 


本 书 特 别 适 合作 为 数学 建 模 竞赛 的 培训 教材 或 参考 用 书 , 也 可 作为 大 学 “数学 实验 ”和 “数学 
建 模 " 以 及 “数据 挖掘" 课程 的 参考 用 书 , 还 可 作为 广大 科研 人 员 、 学 者 ` 工 程 技术 人 员 的 参考 用 书 。 
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作者 简介 


卓 金 武 :硕士 ,2003 年 与 2004 年 获 全 国 大 学 生 数学 建 模 竞 赛 一 等 奖 ,2007 年 获 全 国 研究 
生 数 学 建 模 竞赛 一 等 奖 ,2006 年 获 全 国 研究 生 数学 建 模 竞赛 二 等 奖 ,2004 年 与 2005 年 获 美国 
大 学 生 数 学 建 模 竞赛 二 等 奖 ; 中 国 矿业 大 学 数学 建 模 协会 创始 人 之 一 ,并 担任 第 一 届 数 学 建 模 
协会 执行 主席 ,策划 并 组 织 了 首届 苏 北 高 校 数 学 建 模 联赛 ;多 次 指导 学 生 在 全 国 竞赛 中 获奖 。 


魏 永 生 : 在 读 博 士 , 中 国 工业 与 应 用 数学 学 会 会 员 , 获 国 际 一 等 奖 一 项 ,国际 二 等 奖 两 项 ， 
国家 一 等 奖 两 项 ,二 等 奖 三 项 ,全 国 夏令 营 优 秀 论文 奖 一 项 。 曾 担任 中 国 矿业 大 学 数学 建 模 协 
会 副 主席 , 苏 北 数学 建 模 联 赛 组 委 会 秘书 长 ,第 二 届 苏 北 数学 建 模 联赛 组 织 者 之 一 。 目 前 , 担 
任 北京 交通 大 学 大 学 生 数 学 建 模 竞赛 指导 老师 ,同时 创立 了 北京 诺 亚 数学 建 模 科 技 有 限 公司 
(国内 首 家 提供 数学 建 模 咨 询 服务 的 公司 ,网 址 http://www. noahmodel. com) 。 


秦 健 :硕士 ,中 国 工业 与 应 用 数学 学 会 会 员 , 获 国际 二 等 奖 一 项 ,全 国 大 学 生 数学 建 模 竞赛 
二 等 奖 一 项 ,全 国 研究 生 数学 建 模 竞赛 一 等 奖 一 项 , 曾 担任 中 国 矿业 大 学 数学 建 模 协会 主席 ， 
第 一 届 和 第 二 届 苏 北 数学 建 模 联 赛 组 织 者 之 一 。 现 为 徐州 建筑 职业 技术 学 院 数学 建 模 教练， 
所 指导 的 学 生 在 国际 和 全 国 大 学 生 数 学 建 模 竞赛 中 获得 过 一 等 奖 .二 等 奖 若干 项 。 另 外 ,还 曾 
参与 苏 北 数学 建 模 联赛 的 命题 工作 ,并 多 次 担任 苏 北 数学 建 模 联 赛 评委 。 


李 必 文 : 学 士 , 毕 业 于 安徽 工程 大 学 纺织 服装 学 院 ,2005 年 一 2007 年 先后 在 数学 建 模 大 赛 
中 荣获 安徽 赛区 二 等 奖 、 一 等 奖 和 全 国 一 等 奖 ,2007 年 获 全国 大 学 生 书 画 艺术 大 奖 赛 三 等 奖 。 
曾 任 安徽 工程 大 学 数学 建 模 协会 5200 生 年 安徽 省 十 佳 社团 ) 副 主席 。 自 2008 年 3 月 ,先后 担 
任 华 硕 电脑 (苏州 ) 集 团 PQA 王 程 师 和 南京 匹 尔 害 电 子 信息 技术 有 限 公司 数据 挖掘 工程 师 。 


序 


了 


全 国 大 学 生 数学 建 模 竞 赛 开 展 近 20 年 了 。 在 全 国 竞赛 的 推动 下 ,数学 建 模 教 学 和 竞赛 活 
动 已 经 遍及 国内 大 多 数 高 等 院 校 ,各 种 类 型 的 数学 建 模 教材 .参考 书 、 辅 导 材料 层出不穷 。 这 
些 出 版 物 基本 上 都 出 自 具 有 相当 教学 经 验 的 教师 之 手 。 

与 这 些 书 的 作者 不 同 ,我 应 邀 作 序 并 向 读者 推荐 的 这 本 书 的 编者 是 一 群 <80 后 ”。 他 们 不 
仅 多 次 自己 参赛 并 指导 学 生 参 赛 ,而 且 都 曾 是 学 校 建 模 协 会 、 建 模 联赛 等 活动 的 组 织 者 ,他 们 
通过 参与 竞赛 活动 自身 受益 ,对 数学 建 模 产生 了 很 大 兴趣 ,坚持 不 懈 地 从 事 建 模 的 研究 .教学 
和 应 用 、 推 广 工作 。 他 们 的 这 些 经 历 使 得 这 本 书 有 许多 与 众 不 同 之 处 。 

我 相信 ,这 是 这 群 “80 后 " 们 编写 并 正式 出 版 的 第 一 本 书 , 可 能 会 存在 若干 不 足 。 但 是 , 编 
者 丰富 的 建 模 活动 经 历 反映 到 书 中 ,使 得 一 些 数学 知识 的 应 用 和 处 理 , 一 些 赛 题 的 分 析 和 编 
程 ,特别 是 附件 中 关于 参赛 经 验 的 论述 ,都 具有 鲜明 的 特色 ,值得 准备 参加 竞赛 的 读者 及 数学 
建 模 爱好 者 一 读 。 


姜 启 源 
2010 年 4 月 
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CUMCM( 全 国 大 学 生 数学 建 模 竞赛 ) 已 成 为 大 学 生 参 赛 规模 最 大 的 科技 竞赛 盛事 。 而 
MATLAB 作为 当前 最 优秀 的 数学 软件 ,已 位 列 全 球 最 流行 软件 排行 榜 的 前 20 名 。 在 CUM- 
CM 中 能 否 取 得 好 成 绩 不 仅 取决 于 模型 的 精妙 和 合理 ,还 取决 于 模型 的 求解 。MATLAB 在 模 
型 的 求解 方面 占有 不 可 或 缺 的 位 置 。 目 前 ,CUMCM 设立 了 “MATLAB 创新 奖 ”, 也 说 明 
了 MATLAB 在 数学 建 模 中 的 重要 地 位 。 

MATLAB 入门 很 容易 ,但 要 想 精通 MATLAB 的 所 有 功能 ,几乎 是 不 可 能 的 ,所 以 针对 
自己 的 需要 选择 性 地 学 习 一 些 MATLAB 知识 是 一 种 事半功倍 的 方法 。 目 前 ,关于 MATLAB 
的 书籍 很 多 ,关于 数学 建 模 的 书籍 也 不 少 ,但 将 MATLAB 与 数学 建 模 真 正 结合 在 一 起 的 书籍 
却 不 多 。 本 书 将 从 数学 建 模 的 角度 介绍 MATLAEB 的 应 用 。 

我 们 从 最 初 的 学 生 时 代 自 己 参加 CUMCM 和 MCM( 美 国 大 学 生 数 学 建 模 竞 赛 ) 到 现在 
指导 学 生 参 加 CUMCM 和 MCM, 一 路 走 来 ,有 很 多 体会 。 希 望 通过 这 本 书 , 将 我 们 对 数学 建 
模 和 MATLAB 的 学 习 方法 与 体会 总 结 出 来 ,让 读者 快速 掌握 参加 数学 建 模 竞 赛 应 该 掌握 的 
数学 建 模 方法 和 必 备 的 MATLAB 编程 技术 ,以 提高 竞赛 成 绩 。 在 内 容 的 编排 上 ,也 是 尽量 做 
到 让 初学 者 能 够 大 门 , 让 经 验 丰富 者 能 耳目 一 新 ,思路 胃 然 。 书 中 融入 了 我 们 许多 参加 CUMCM 
和 MCM 的 真实 体会 以 及 多 年 一 线 指导 学 生 参 加 竞赛 的 经 验 ,使 本 书 具 有 实用 、 创 新 和 鲜 活 的 
特征 。 

全 书 分 上 下 两 篇 讲述 MATLAB 在 数学 建 模 中 的 应 用 。 上 篇 中 第 1 一 3 章 讲 的 是 数学 建 
模 中 常规 方法 的 MATLAB 实现 ,读者 看 完 这 部 分 后 ,就 可 以 用 MATLAB 解决 数学 建 模 中 的 
一 些 基 础 问题 了 ,如 数据 处 理 、. 程 序 绘图 ,求解 目标 规划 模型 .数据 量 较 少时 的 预测 等 。 随 着 数 
学 建 模 竞赛 的 发 展 , 题 目 对 参赛 者 创造 性 解决 问题 的 要 求 越 来 越 高 。 也 就 是 说 ,参赛 者 需要 有 
自主 开发 程序 的 能 力 , 才 有 可 能 在 比赛 中 取得 好 成 绩 * 所 以 就 要 学 习 更 为 复杂 的 编程 方法 和 技 
巧 。 于 是 在 上 简 剩 下 的 简 帆 安排 了 数学 建 模 高 级 方法 的 MATLAB 实现 ,包括 遗传 算法 ,粒子 
群 算法 、 模 拟 退 火 算法 .神经 网 络 ,小波 分 析 ` 动 态 模拟 、 数 值 模 拟 等 。 通 过 这 些 内 容 的 学 习 , 读 
者 基本 上 就 可 以 求解 一 些 复杂 的 问题 了 。 现 在 的 建 模 竞 赛 经 常会 出 现 的 问题 不 是 模型 建 不 出 
来 ,而 是 模型 求解 不 出 来 ,所 以 高 水 平 的 MATLAB 编程 和 对 常用 方法 的 熟悉 将 对 参加 数学 建 
模 竞 赛 相当 有 帮助 。 下 篇 介绍 了 数学 建 模 竞赛 中 典型 真题 的 MATLAB 求解 方法 及 过 程 , 通 
过 对 一 些 真题 案例 的 解析 ,加 强 读 者 对 MATLAB 的 驾驭 能 力 , 巩 固 读 者 对 MATLAB 的 操作 
能 力 ,同时 增加 读者 对 数学 建 模 问题 的 敏感 度 。 

纵 观 全 书 , 可 发 现 该 书 的 特点 相当 鲜明 ,主要 表现 在 如 下 几 个 方面 

(1) 结构 编排 合理 。 在 讲解 方法 时 ,由 浅 和 人 深 , 循 序 渐进 ,让 初学 者 知道 人 门 的 切入 点 ,让 
建 模 和 编程 高 手 又 有 借鉴 的 空间 。 书 中 有 很 多 小 程序 可 供 初 学 者 套用 及 学 习 , 又 不 乏 实现 功 
能 复杂 的 篇 幅 较 长 的 程序 ,这 部 分 主要 面向 学 有 余力 的 读者 。 

(2) 既 突 显 数学 建 模 思想 又 包含 了 完整 的 MATLAB 程序 。 目 前 , 现 有 的 数学 建 模 书 籍 
完整 附带 程序 的 较 少 ,而 本 书 给 出 了 大 量 实 用 、 完 整 . 详 细 的 源 代 码 , 这 对 于 初学 数学 建 模 和 
MATLAB 的 读者 来 说 ,具有 非常 大 的 参考 价值 。 

(3) 部 分 新 型 方法 是 同类 书籍 中 比较 少见 的 。 本 书 中 的 动态 仿真 和 小 波 分 析 是 目前 其 他 
建 模 书籍 中 比较 少见 的 。 动 态 仿真 ,美国 大 学 生 刚 开始 尝试 在 建 模 比赛 中 使 用 ,在 国内 尚 属 空 











白 ; 小 波 分 析 虽 有 些 难 懂 ,但 确实 是 一 种 有 力 的 工具 ,在 本 书 中 也 对 其 进行 了 深入 浅 出 的 介绍 。 

《4) 本 书 涉及 知识 面 宽广 ,以 建 模 为 中 心 ,辐射 计算 机 、 生 物 学 、 环 境 科 学 ,心理 学 .古典 文 
学 .影视 .服装 设计 以 及 艺术 学 。 比 如 ,将 污染 物 扩 散 的 四 维 水 质 模型 与 服装 设计 面料 的 透气 
透 湿性 能 联系 起 来 , 找 出 信息 隐藏 的 数字 水 印 技术 与 (孙子 兵法 ) 中 的 “出 其 所 不 趋 , 趋 其 所 不 
意 ” 的 思想 结合 ,等 等 。 书 中 所 涵盖 的 这 些 领 域 虽然 只 是 被 略微 带 过 ,但 从 侧面 已 折射 出 数学 
建 模 是 广泛 服务 于 社会 各 个 领域 的 。 在 现代 社会 , 某 学 科 单打 独 斗 的 时 代 已 经 过 去 了 ,本 书 在 
无 形 之 中 已 经 树立 了 一 个 意识 :知识 之 间 是 相通 的 。 

〈5) 附件 中 关于 参加 建 模 竞赛 的 经 验 体会 也 在 同类 书籍 中 比较 少见 ,也 都 是 我 们 参加 数 
学 建 模 竞赛 的 真实 体会 。 我 们 站 在 一 个 成 功 参赛 者 的 角度 与 读者 进行 交流 ,更 容易 让 读者 感 
同 身受 。 比 如 , 书 中 介绍 了 数学 建 模 中 的 项 目 管理 和 时 间 管 理 ,这 些 都 是 容易 被 参赛 队员 忽视 
却 对 竞赛 成 绩 影响 相当 大 的 因素 。 再 比如 , 书 中 介绍 了 目标 建 模 方法 ,这 是 一 种 非常 实用 的 数 
学 建 模 方法 。 有 了 这 些 经 验 和 实用 方法 的 指导 ,读者 会 有 的 放 矢 地 去 准备 数学 建 模 竞 赛 , 同 时 
也 能 增强 读者 参加 数学 建 模 竞 赛 的 信心 。 

《6) 书 中 内 容 以 数学 建 模 竞赛 为 主线 ,同时 有 意 将 数学 建 模 向 实际 的 生产 力 方向 推广 , 真 
正体 现 了 数学 建 模 的 魅力 。 这 具体 体现 在 第 2 章 规划 问题 所 涉及 的 投资 问题 ,第 10 章 彩 票 问 
题 以 及 第 14 章 出 版 社 资源 配置 问题 等 内 容 的 编排 上 。 建 模 如 果 只 是 单纯 的 论文 ,不 能 产生 社 
会 价值 ,那么 建 模 永 远 只 会 停留 在 学 术 层面 ,不 能 创造 社会 价值 就 等 于 没有 价值 。 所 以 纵 观 通 
篇 书稿 ,本 书 都 在 有 意 引导 读者 思考 如 何 让 建 模 产生 更 实际 的 社会 价值 。 

《7) 书 中 原创 知识 占据 绝 大 部 分 版 面 ,让 读者 少 了 儿 分 似曾相识 或 模型 雷同 的 感受 。 常 
规 方法 \ 高 级 方法 和 真题 案例 的 结构 部 署 也 让 本 书 独 树 一 帜 ,让 读者 在 学 习 MATLAB 和 数学 
建 模 过 程 中 有 一 个 循序 渐进 的 过 程 ;使 读者 在 短 时 间 内 成 为 一 个 MATLEAB 程序 开发 高 手 , 同 
时 成 为 一 个 数学 建 模 高 手 。 

《8) 内 容 详实 ,文字 活泼 ,字里行间 透射 灵动 鲜 活 的 魅力 。 在 本 书 编写 过 程 中 ,在 保证 描 
述 精准 的 前 提 下 ,我 们 据 弃 了 那些 刻板 ,索然 无 味 的 文字 ,让 文字 也 有 了 活力 。 
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杭 所 退火 | 生生 于 全 和 和 生机 划 | 态 ,下 载 ASA 程序 包 并 深度 解读 ,让 初学 
人 者 能 入 门 , 让 热 知 该 算法 的 读者 思路 项 然 











开朗 








续 表 





方法 名 称 


在 数学 建 模 中 的 应 用 场合 


特色 内 容 





人 工 神 
经 网 络 


数学 建 模 中 一 切 聚 类 .评价 及 预 
测 模式 的 问题 


@ 在 理论 上 并 列 解析 BP 网 络 和 RBF 
网 络 ; 

加 从 全 新 角度 阐述 神经 网 络 的 基本 原理 ， 
避免 与 其 他 同类 书籍 雷同 ,偏重 原创 和 
细致 ; 

图 设计 四 层 径 向 基 小 波 神经 网 络 结构 ; 
图 同时 用 MATLAB 神经 网 络 工具 箱 和 开 
发 MATLAB 源 程序 对 公路 运 量 问题 进行 
预测 








小 波 分 析 


四 海量 数据 趋势 挖掘; 
四 组 建 小 波 神经 网 络 


@ 多 角度 切 人 小 波 神经 网 络 ; 
加 引入 融合 数学 建 模 思想 和 小 波 分 解 思想 
的 数字 水 印 技术 





计算 机 虚拟 








四 模型 难以 求 出 解析 解 ; 
加 动态 (动画 ) 展 现 数学 模型 
图 复杂 非 线 性 规划 问题 粗 略 求解 





〇 有 趣 解说 求解 的 蒙特 卡 罗 模 拟 # 

回 动态 仿真 表现 力 强 ?但 是 当前 具有 动态 
仿真 的 数学 建 模 书籍 还 不 多 见 ， 

图 用 四 种 蒙 氏 模拟 法 解决 非 线性 规划 
轿 电 影 (功夫 熊猫 ) 武 打 片段 电影 仿真 ， 
回 完整 地 组 建 .求解 四 维 水 质 模型 ,并 可 用 
计算 机 模拟 连续 的 瞬间 水 质 情形 








第 工 章 
| 数据 建 模 常规 方法 的 MATLAB 实现 


数据 处 理 是 数学 建 模 的 基础 ,通常 遇 到 的 问题 是 对 采集 到 的 数据 进行 处 理 和 分 析 , 从 而 得 
到 这 些 数据 所 反映 的 信息 。 从 数学 建 模 的 角度 ,将 数据 反映 出 来 的 信息 转化 成 数学 表达 式 是 
建 模 的 基础 ,所 以 通常 对 数据 的 处 理 就 是 趋势 分 析 和 将 数据 转化 为 函数 表达 式 , 这 里 通常 用 到 
的 是 数据 的 插值 和 曲线 拟 合 。 应 该 说 ,MATLAEB 的 应 用 使 得 这 些 方法 很 容易 实现 ,尤其 当 数 
据 量 较 大 时 ,MATLAB 的 数据 处 理 优 势 就 更 明显 了 。 


1.1 数据 的 读 入 与 读 出 








用 户 用 MATLAB 进行 科学 计算 时 不 可 避免 地 要 用 到 大 量 的 数据 ,而 方便 的 数据 处 理 方 
法 会 让 用 户 更 得 心 应 手 。MATLAB 提供 了 多 种 处 理 数据 的 方法 ,一 种 方法 是 将 数据 输出 , 然 
后 复制 粘贴 到 其 他 软件 中 进行 处 理 , 但 是 这 种 方法 不 方便 ; 勇 一 种 方法 是 与 Excel 和 记事 本 进 
行 数据 交互 ,这 种 方法 方便 了 数据 的 操作 。 


1.1.1 Excel 与 MATLAB 的 数据 交互 


MATLAB 中 的 Excel Link 工具 是 一 个 实现 与 Excel 进行 交互 的 插件 。 通 过 连接 Excel 
和 MATLAB, 用 户 可 以 在 Excel 工作 表 空 间 和 安 编 程 工具 中 使 用 MATLAB 的 数值 计算 和 图 
形 处 理 等 功能 ,不 需要 脱离 Excel 环境 :同时 由 Excel Link 保证 两 个 工作 环境 中 的 数据 交换 
和 同步 更 新 。 

Excel Link 的 安装 :首先 ;要 确保 系统 中 正确 安装 了 Excel 软件 .MATLAB 和 Excel Link 
[ 具 箱 。 启 动 Excel, 选 择 荣 单 “工具 ?项 下 的 “加 载 宏 " 项 ,弹出 一 个 加 载 宏 对 话 框 ,选中 Excel 
Link 项 (一 般 是 Spreadsheet Link EX for use with MATLAB)。 如 果 该 项 不 存在 , 则 通过 浏 
览 目录 ,在 目录 MATLABNtoolbox\exlink 文件 夹 下 找到 excllink. xla 文件 , 单 击 “ 确 定 ”, 选 中 
Excel Link 项 并 确定 后 ,在 Excel 中 多 了 一 个 Excel Link 工具 条 ,如 图 1-1 所 示 。 


j startmat lab Dputmatrix getmatrix evalstring &etfigure Wizard preferemnces 肌 


图 1-1 Excel Link 工具 条 


经 过 以 上 的 设置 后 就 可 以 使 用 Excel Link 了 ,读者 看 到 工具 条 后 会 很 容易 学 会 其 操作 ， 
这 里 不 再 袭 述 ,有 问题 可 参考 Help。 


1.1.2 记事 本 与 MATLAB 的 数据 交互 


在 MATLAB 中 可 以 读 取 存 人 记事 本 中 的 数据 ,也 可 以 将 内 存 中 的 变量 数据 保存 到 记事 
本 中 ,以 便 后 期 操作 。 








1. MATLAB 读 取 记 事 本 中 的 数据 


当 记 事 本 中 记录 的 全 部 都 是 数据 时 ,可 以 用 函数 load, 其 形式 是 : 


load('filename. # 关 关 7) 


其 中 ,文件 扩展 名 可 以 是 任意 记事 本 文件 的 扩展 名 ,运行 此 函数 后 就 会 把 记 寻 


按 和 矩阵 的 形式 放 人 名 为 filename 的 变量 中 。 








本 中 的 数据 


当 记事 本 中 的 数据 结构 变 得 复杂 时 ,例如 表 1- 1 中 的 数据 (第 1 行 除外 ), 就 不 能 用 本 数 


load 了 ,此 时 函数 textread 是 最 优 的 选择 ,其 规则 是 : 


[LA,B,C,…]=textread('filename”，format”,N) 





其 中 ,A,B,C,… 为 每 一 列 数据 将 要 保存 的 变量 名 ;format 为 读 取 格式 ;N 为 读 取 次 数 。 
例如 上 面 的 文件 保存 为 当前 路 径 下 的 t+. txt 文 件 后 ,输入 命令 得 到 : 


>> [name,type,x,Yyanswer] = textread( t.txt'，$S TYPeg%n 多 开 mn 有 5 2) 

















name = 
"Sally 
二 
type = 
x = 
3.1000 
2.5000 
Y = 
45 
20 
answer = 
Yes 
7 
表 1- 工 数 据 
names types z y answer 
Sally Typel 3.1 45 Yes 
Tom Typel 2.5 20 No 
joy Type3 0.51 1 No 
Lily Type2 2.1 80 Uncertain 























在 命令 中 ' 失 s Type%n 5 5%n 5s? 即 为 format 的 内 容 , 最 后 一 个 参数 2 表示 读 取 2 行 


数据 。 


在 MATLAB 中 ,还 有 一 个 函数 fscanf 可 以 读 取 记 寻 


需要 ,请 直接 在 Help 中 寻求 帮助 。 
2. 将 MATLAB 数据 写 入 记事 本 中 





本 中 的 文件 ,其 功能 更 加 强大 ,如 有 


同样 ,如 果 只 需要 保存 数据 ,并 且 以 后 也 只 用 MATLAB 调用 ,可 以 用 save 命令 ,其 语 


法 是 : 
save file objl obj2… 


将 各 变量 objl,obj2… 存 人 文件 名 e 中 (为 . mat 格式 ) ,用户 无 法 直接 查看 文件 中 保存 





的 数据 。 

如 果 想 将 数据 保存 为 平常 的 . txt 文档 ,可 以 由 其 他 软件 进行 读 取 , 函数 fprintf 是 一 个 不 
错 的 选择 。 还 以 上 面 的 数据 为 例 , 要 将 这 些 量 按 格式 保存 在 一 个 tp. txt 文档 中 ,可 以 建立 以 
下 命令 文件 。 


fid = fopen( tp.txt'，wt )5 

fprintf(fid， This is the database of class 1.\n ); 

name = "Sally';y types=1; x=3.1;Y=45; answer = Yes 
fprintf(fid,$S Type$%u 各 上 $u 名 5 \n' ynameytypes,xyyyanswer) 1 
name = "Tom'; types=1; Xx=2.5;Y=20; answer= No 
fprintf(fid, 多 s Typeg%u 要 上 $u 有 5 \n'ynameytypes,xyyyanswer); 
name =' Joy'; types = 3; x= .51; Y=1; answer= No 

fprintf(fid, % s TYpe 和 各 上 名 $S \Nn nameytypes,xyyyanswer); 
name = "Lily'; types=2; x=2.1;Y=80; answer=' Uncertain 
fprintf(fidy 多 S Type%u 和 上 多 和 5 \n ynameytypes,xyyyanswer); 
fclose(fid) 


运行 后 ,就 会 建立 一 个 tp. txt 文档 ,其 中 的 内 容 如 下 ， 


This is the database of class 1 . 
Sally Typel 3.100000 45 Yes 

Tom Typel 2.500000 20 No 

Joy Type3 0.510000 1 No 

Lily TYpe2 2. 100000 80 Uncertain 


现在 对 程序 中 的 语句 进行 说 明 。 第 1 段 是 打开 一 个 文件 ,返回 一 个 指标 fid, 第 2 段 是 直 
接 在 文件 中 写 人 一 名 话 “This is the database of class 1. ”, 然 后 就 是 每 赋值 一 次 ,就 向 文件 中 
写 人 一 次 : 


fprintt(fjd, 和 sypeyu wyf% 了 可 5%5\nTynameitypes,, 册 answer)y 
其 格式 为 
fprintf(id， format ,argl,arg2，…)5 


其 中 ,%s 表示 字符 串 (string); %f 表示 浮 点 数 (float) ; %u 表示 十 进 制 数 。 

另外 还 有 很 多 格式 ,用 时 请 参阅 Help 中 的 fprint 以 及 其 下 的 formatting strings。 

在 这 里 还 要 注意 ,不 能 用 郴 数 textread 读 取 后 再 用 郴 数 fprintf 进行 写 操作 。 这 是 因为 本 
数 textread 读 取 数据 后 保存 为 cell 结构 ,而 函数 fprintf 不 能 进行 cell 结构 数据 的 写 操作 。 


1.2 数据 拟 合 方法 


曲线 拟 合 也 称 为 曲线 逼近 , 它 与 揪 值 函数 有 些 区 别 , 只 要 求 拟 合 的 曲线 能 合理 地 反映 数据 
的 基本 趋势 ,而 并 不 要 求 曲线 一 定 通 过 数据 点 。 曲 线 拟 合 有 几 种 不 同 的 判别 准则 ,如 使 偏差 的 
绝对 值 之 和 最 小 、 使 偏差 的 最 大 绝对 值 最 小 和 使 偏差 的 平方 和 最 小 ( 即 最 小 二 乘法 ) 。 常 用 的 
方法 是 最 后 一 种 。 


1.2.1 多 项 式 拟 合 


将 数据 点 按 多 项 式 的 形式 进行 拟 合 ,使 用 最 小 二 乘法 可 以 确定 多 项 式 的 系数 。 多 项 式 拟 
合 有 指令 语句 和 图 形 窗 口 两 种 方法 。 


1. 多 项 式 拟 合 指令 
polyfit(X,Y,N) :多 项 式 拟 合 ,返回 降 索 排 列 的 多 项 式 系数 。 
polyval(P,xi) :计算 多 项 式 的 值 。 
其 中 ,X,Y 是 数据 点 的 值 ;N 是 拟 合 的 最 高 次 宕 ;P 是 返回 的 多 项 式 系数 ;xi 是 要 求 的 点 
的 横 坐标 。 
例如 要 对 表 1 - 2 中 的 数据 进行 多 项 式 拟 合 。 
表 1-2 实例 数据 








1 2 3 4 5 6 7 8 9 






































y 9 7 6 3 | 2 5 7 20 





拟 合 的 命令 如 下 : 


x=[123456789]; 
Y=[9763 -125720]; 
P= polyfitCxyy,3)5 

xi= 0:.2:10 

吐 = polyval(P,xi)4 
plot(xiyyiyx,yy， rw ) 5 


得 到 的 拟 合 曲线 与 原始 数据 的 对 照 如 图 1-2 所 示 。 








图 1-2 拟 合 曲线 与 原始 数据 的 对 照 


2. 图 形 窗口 的 多 项 式 拟 合 
在 图 形 窗 口中 可 以 用 菜单 方式 对 数据 进行 简单 快速 高效 的 拟 合 。 以 表 1 - 2 中 的 数据 
为 例 , 先 画 出 数据 点 : 





>>x=[123456789]iy=[9763 -125720];plot(x,yyrx 7)i; 





然后 在 图 形 窗口 单 击 Tools-~Basic Fitting, 打 开 对 话 框 ,并 按 图 1 - 3(a) 进 行 操作 ,就 可 以 得 


到 图 1- 3(b) 的 效果 ,其 中 分 别 用 线性 .二 阶 、 三 阶 对 数据 进行 了 多 项 式 拟 合 ,下 


示 残 差 ,可 以 看 出 ,进行 三 阶 多 项 式 拟 合 的 效果 是 最 好 的 ,与 前 述 结果 一 样 。 


CTTTTEE 和 | 


Select data: EECED 全 天 国 ~ 
回 Center and scale x data 
Plot fits 
Check to _ display fits on 人 ur 
spline interpolant 
Shape-preserving interp... 
] linear 
] quadrat ic 
cubie 
外 h degree polynomial 
5th degree polynomial 
6th degree Polynomial 
回 Show equat ions 


Sisnificant disgits:12 “ 


回 Plot residuals 


Bar plot > 
Subplot ~ 
口 Show norm of residuals 


Help Close 





图 1-3 图 形 窗口 的 多 项 式 拟 合 


1.2.2 指定 函数 拟 合 
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在 MATLAB 中 也 可 以 用 用 户 自 定义 的 函数 进行 拟 合 。 通 过 下 面 的 例子 可 以 了 解 用 指定 


函数 进行 数据 拟 合 的 基本 方法 。 


在 某 次 阻尼 振 划 实验 中 测 得 如 表 1- 3 所 列 的 18 组 数据 点 。 
表 1-3 数据 点 的 坐标 












































灾 0 0.4 1.2 2 2.8 3.6 4.4 5.2 6 
3 1 0.85 | 0%.29 | 一 0.27 | 一 0.53 | 一 0.4 | 一 0.12| .17 | %328 
7.2 8 9.2 104 | 1.6 | 124 | 13.6 | 14.4 15 
> 0.15 | 一 0.03 | 一 0.15 | 一 0.071| 0.059 | 0%.08 | o.032 | 一 0.015| 一 0.02 








在 MATLAB 中 可 以 用 如 下 命令 画 出 上 述 点 的 散 点 图 ( 





见 图 1- 4): 





0.4;51.2;2;2.8;3.6;34.4;5.2;6;7.238;9.2;10.4;11.6;12.4;13.6;14.4315]; 

Y= [1;0.85;0.29; - 0.27; -0.53; -0.4; 一 0.12;0.17;0.28;0.15; -0.03; -0.15; 一 0.071;0.0593 
0.08;0.032; - 0.015; - 0.02]; 

plot(Cxyy， 工 关 )5 

















5 10 1 


图 1-4 数据 散 点 图 


知道 其 对 应 的 函数 形式 为 FlO) =acosCkt)e"… , 则 可 用 MATLAB 进行 拟 合 。 编 写 如 下 的 
M 文 件 ; 


Syms 七 
x= [0;0.4;1.2;212.813.614.415.216;7.2384972310.4;11.6112.4113.6114.4;15] 
Y= [1;0.85;0.29;- 0.27;- 0.53; 一 久生 一 0.12;0.1740.28;0.15; - 0.03; -0.15; - 0.071;0.059; 
0.08;0.032; -0.015; -0.02]; 
= fittype( ax cos(Kx t) x exp(W 汪 it ，independent ，t'，coefficients ,ayk yw)); 
efua=Eit(x,Wt) #% 显示 拟 合 函 数 
邓 E 0 20 
对 = cfun(xi)y 
plot(x,yy，rw yxiyyiyb=); 


运行 此 程序 ,在 Command Window 可 得 到 如 下 的 运行 结果 (图 形 结果 如 图 1- 5 所 示 ): 


Warning: Start point not provided，choosing random start point. 
> Infit > handlewarn at 714 
In fit at 328 
cfun = 
General model: 
cfun(t) = axcos(kxt)x expCwxt) 
Coefficients (with 95#% confidence bounds) : 


a = 0.9987 《0.9836，1.014) 
:= 1.001 《0.9958， 1.006) 
= -0.2066 〈-0.2131，- 0.2002) 


程序 中 ,函数 fittype 是 自 定 义 拟 合 函数 ;cfun 一 fit(x,y, 旭 是 根据 自 定 义 的 拟 合 函数 了 来 
拟 合 数据 z,y。 注 意 :此 处 数据 必须 为 列 向 量 的 形式 。 画 数 fittype 和 枯 数 fit 的 用 法 与 参数 


规则 可 参考 Help。 

从 结果 可 以 看 出 , 拟 合 的 曲线 为 F(z) 一 0. 998 7cos(1. 001t)e “2 。 拟 合 曲线 给 出 了 数 
据 的 大 致 趋势 ,结果 很 好 ,并 给 出 了 各 参数 的 置信 区 间 。 

注意 :Command Window 中 给 出 的 Warning, 是 由 aze 3 个 参数 的 初始 值 未 给 出 导致 
的 ,因此 如 果 拟 合 结果 不 理想 ,可 以 多 运行 几 次 。 


1.2.3 曲线 拟 合 工具 箱 


MATLAB 提供 了 曲线 拟 合 工具 箱 , 其 功能 很 强大 , 使 用 也 很 方便 。 依 次 单 击 Start 一 
Toolboxes-~Curve Fitting-Curve Fitting Tool(Ccftool) 就 可 以 打开 曲线 拟 合 工具 箱 ,或 者 在 
Command Window 中 直接 输入 cftool 命令 打开 。 打 开 后 的 界面 如 图 1-6 所 示 。 
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图 1-5 程序 运行 结果 图 1-6 曲线 拟 合 工具 箱 界面 


界面 中 5 个 按钮 的 功能 分 别 如 下 : 

Data: 输 出 查看 和 平滑 数据 : 

Ritting: 拟 合 数据 ,比较 拟 合 山 线 和 数据 焦 。 

Exclude: 可 从 拟 侣 曲线 中 排除 特殊 的 数据 点 ; 

Plotting : 选 定 区 间 后 * 单 击 按钮 可 选择 原始 数据 和 拟 合 的 数据 作 图 。 

Analysis: 对 拟 合 进 行 满意 度 .偏差 等 分 析 。 

曲线 拟 合 工具 箱 中 包含 了 各 种 常用 的 数据 拟 合 方法 ,可 以 对 各 种 函数 进行 拟 合 ,具体 的 操 
作 方 法 请 参考 Help。 


1.3 数据 拟 合 应 用 实例 


1.3.1 人 口 预测 模型 


表 1-4 是 某 地 区 1971 一 2000 年 的 人 口 数据 , 试 给 出 该 地 区 人 口 增长 的 数学 模型 。 根 据 
表 中 的 数据 做 出 如 图 1- 7 所 示 的 散 点 图 。 

从 图 1- 7 可 以 看 出 ,人 口 随时 间 的 变化 呈 非 线性 过 程 , 而 且 存在 一 个 与 横 坐 标 轴 平 行 的 
渐 近 线 , 故 可 以 用 Logistic 曲线 模型 进行 拟 合 。 





表 1-4 某 地 区 人 口 变化 数据 












































年 份 | 时 间 变 基 ! 一 (! 年份 一 1970) | 人 口 >/ 人 |‖ 年 份 | 时 间 变 基 : 一 (! 年 份 一 1970)》 | 人 口 >/ 人 
1971 1 33815 1986 16 34520 
1972 2 33981 1987 17 34507 
1973 3 34004 1988 18 34509 
1974 4 34165 1989 19 34521 
1975 5 34212 1990 20 34513 
1976 6 34327 1991 21 34515 
1977 7 34344 1992 22 34517 
1978 8 34458 1993 23 34519 
1979 日 34498 1994 24 34519 
1980 10 34476 1995 25 34521 
1981 11 34483 1996 26 34521 
1982 12 34488 1997 2 34523 
1983 13 34513 1998 28 34525 
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图 1-7 某 地 区 人 口 随时 间 变 化 的 散 点 图 
因为 Logistic 曲线 模型 的 基本 形式 为 








一。 
>? 4 十/ 

所 以 ,只 要 令 ,就 可 以 将 其 转化 为 直线 模型 : 
一 aa 十 ia 


下 面 分 别 用 MATLAB 进行 回归 分 析 拟 合计 算 。 首 先 ,编写 MATLAB 回归 拟 合 程序 ， 
见 P1 一 1。 





程序 编号 | P1-1 | 文件 名 称 | Population. m | 说 明 人 口 Logistic 曲线 模拟 
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clear 
cle 
#% 读 人 人 口 数 据 (1971- 2000) 
Y= [33815 ”33981 “34004 ”34165 ”34212 ”34327 ”34344 。 34458 。 34498 。 34476 
34483 ”34488 。 34513 ”34497 ”34511 ”34520 ”34507 ”34509 ”34521 。 34513 
34515 ”34517 ”34519 ”34519 ”34521 。” 34521 ”34523 ”34525 。” 34525 。 34527] 
% 读 人 时 间 变量 数据 (t= 年 份 -1970) 
mT= [1 当 3 竺 对 6 旦 8 9 10 11 12 13 14 15 
7 18 19 20 21 22 23 24 25 26 2 28 29 30] 
# 线性 化 处 理 
fort = 1:30， 
x(t) = exp( 一 tb 
YCt) = 1/YCt)5 
end 
% 计算 ,并 输出 回归 系数 B 
ec= zeros(30,1) +15 
X= [cx]; 
B= invCX' we X) a 和 于 
for i=1:30， 
% 计算 回归 拟 合 值 
zGD =B(Gl,1) +B(2,1) wx(iD5 
% 计算 离 差 
Si) =Y(Gi) - sum(y)/303 
% 计算 误 郑 
wGD =z(G-YCGDi 
end 
#% 计算 离 差 平方 和 S 
S=syws4 
#% 回归 误差 平方 和 0 
Q5Wx WE 
多 计算 回归 平方 生 
U=s-0; 
s% 计算 ,并 输出 了 检验 值 
F=28x 0U/9 
s#% 计算 非 线性 回归 模型 的 拟 合 值 
for j=1:30， 
Y(G) = 1/(BG1,1) + B(2,1) * exp(-j)); 
end 
s* 输出 非 线性 回归 模型 的 拟 合 曲 线 (Logisic 出 线 ) 
plot(T,Y) 








运行 该 程序 后 ,输出 ; 


(1) 回归 系数 日 及 F 检验 值 : 





1.0e-004 * 
0.2902 
0. 0182 


47. 8774 


(2) Logistic 模型 拟 合 曲线 ( 见 图 1- 8)。 
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图 1-8 Logistic 模型 拟 合 曲线 


1.3.2 薄膜 渗透 率 的 测定 


1， 问 题 的 提出 

某 种 医用 薄膜 有 允许 一 种 物质 的 分 子 穿 透 它 (从 高 波 
度 的 溶液 向 低 浓度 的 溶液 扩散 ) 的 功能 ,在 试制 时 需 测定 薄 
膜 被 这 种 分 子 穿 透 的 能 力 。 测 定 方法 如 下 :用 面积 为 $ 的 
薄膜 将 容器 分 成 体积 分 别 为 VA Vs 的 两 部 分 ( 见 图 1- 9)， 
在 两 部 分 中 分 别 注 满 该 物质 的 两 种 不 同 浓度 的 溶液 。 此 
时 ,该 物质 分 子 就 会 穿 过 薄膜 从 高 浓度 溶液 向 低 浓度 溶液 
扩散 。 已 知 通过 单位 面积 薄膜 分 子 扩散 的 速度 与 膜 两 侧 溶 
液 的 浓度 差 成 正比 ,比例 系数 民 表 征 了 薄膜 被 该 物质 分 子 图 1-9 圆柱 体 容器 被 薄膜 





穿 透 的 能 力 , 称 为 渗透 率 二 定时 测量 容器 中 薄膜 某 一 侧 的 敌 面 阻 忆 
溶液 浓度 值 ,可 以 确定 K 的 数值 ,试用 数学 建 模 的 方法 解决 K 值 的 求解 问题 。 

2. 数学 模型 的 建立 

(1) 假 设 


薄膜 两 侧 的 溶液 始终 是 均匀 的 , 即 在 任何 时 刻 膜 两 侧 的 每 一 处 洲 液 的 浓度 都 是 相 
同 的 。 

Q@ 当 两 侧 浓度 不 一 致 时 ,物质 的 分 子 穿 透 薄膜 总 是 从 高 浓度 溶液 向 低 浓度 溶液 扩散 。 

@ 通过 单位 面积 膜 分 子 扩散 的 速度 与 膜 两 侧 溶液 的 浓度 差 成 正比 。 

@@ 薄膜 是 双向 同性 的 , 即 物质 从 膜 的 任何 一 侧 向 另 一 侧 渗 透 的 性 能 是 相同 的 。 

(2) 符号 说 明 

@ CA(D .Cs(G) 表 示 上 时 刻 膜 两 侧 溶 液 的 浓度 。 

四 aa ,as 表示 初始 时 刻 膜 两 侧 溶 液 的 浓度 〈 单 位 :mg/cnm ) 。 

加 K 表示 渗透 率 。 

@ WA Vs 表示 由 薄膜 阻隔 的 容器 两 侧 的 体积 。 





(3) 模型 的 建立 

考察 时 段 [:,z 十 Ai] 薄 膜 两 侧 容 器 中 该 物质 质量 的 变化 。 以 容器 A 侧 为 例 , 在 该 时 段 物 质 
质量 增加 量 为 

VaACAG 十 At) 一 VACACtD) 
另 一 方面 由 渗透 率 的 定义 可 知 , 从 B 侧 渗透 至 A 侧 的 该 物质 的 质量 为 
SK (Cs 一 CA)At 
由 质量 守恒 定律 ,两 者 应 该 相等 ,于 是 有 
VACACG 十 A) 一 VACA(C) 一 SK(Csa 一 CA)A 

两 边 同 除 A, 令 At 一 0, 整 理 得 





只 = 盖 cw G-D 
且 注 意 到 整个 容器 的 溶液 中 含有 该 物质 的 质量 应 该 不 变 , 即 有 下 式 成 立 ， 
VACACD) 十 VaCaG) 一 VAaA 十 Vaas CG1-2a) 
CaGD = ax + 更 w 一 草 Ca) (1-2b) 
了 本 
将 式 (1- 2b) 代 入 式 (1- 1) 得 
王国， TS 
是 +sKk{( 寺 + 直径 = 区 作 攻 站 六 
再 利用 初始 条 件 Ca(0) = 一 oz 
二 4AVA 十 aaVs | VA(GB 一 aA) - 呈 ( 上 + 吉 )， 
解 得 GD 二 SR 秆 这 (R+ 古 ) 


至 此 ,问题 归结 为 利用 Cs 在 时 刻 5 的 测量 数据 C, 9 =1,2,…，,N) 来 辨识 参数 K 和 aA、 
aa ,对 应 的 数学 模型 变 为 求 本 数 ， 


minE(KYa 虽 ,as) 一 > (Cs 如) 一 CD) 


产 
全 和] 【加 作 一 和 
二 这: 到 下 了 





问题 转化 为 求 本 数 
ECKvanraa) 一 局 |oHocs st 一 C] 


的 最 小 值 点 (Ka 0) 。 
3. 求解 参数 
例如 , 设 VA=Va 一 1000em' ,S 一 10cm ,对 容器 了 部 分 溶液 浓度 的 测试 结果 如 表 1- 5 所 列 。 
表 1-5 容器 B 部 分 溶液 测试 浓度 





三 100 200 300 400 500 600 700 800 900 1000 









































CiV/(mg em-3)| 454 | 499 | 5.35 | 5.65 | 5.90 | 6.10 | 6.26 | 6.39 | 6.50 | 6.59 





此 时 极 小 化 的 函数 为 
ECK,aaAyas) 一 六 [5 十 pe ok 一 C 〇 了 
下 面 用 MATLAB 软件 进行 计算 。 


(1) 编写 M 文件 (curvefun. m) 


function 上 = curvefun(xvtdata) 
三 = x(1) +Xx(2) * exp( 一 0.02x*x(3) * tdata); 当 其 中 ,x(1) = ax(2) = b,x(3) = k 


(2) 编写 程序 (testl. m) 


tdata = linspace(100,1000,10)3 

cdata = le- 05. * [454 499 535 565 590 610 626 639 650 659]; 
x0 = [0.2,0.05,0.05]; 

x = 1sqcurvef 让 (〔 curvefun ,x0,tdataycdata) 

上 = curvefun(xtdata) 

plot(tdatavcdata， ovtdata;E， 工 一 ) 


(3) 输出 结果 
x = 
0.0063 ”一 0.0034 0.2542 
即 表 示 k = 0.2542， a = 0.0063， b = 一 0.0034 
上 = 


0.0043 ”0.0051 ”0.0056 0.0059 “0.0061 0.0062 ”0.0062 ”0.0063 
0.0063 0.0063 


曲线 的 拟 合 结果 如 图 1= 10 所 示 , 进 一 步 可 求 得 本 一 0.004 mg/cms5aa 一 0.01 mgjem' 。 
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图 1- 10 模型 拟 合 曲线 与 溶液 实际 测试 浓度 


1.4 数据 的 可 视 化 


1.4.1 地 形 地 和 狐 图 形 的 绘制 


表 1-6 所 列 为 某 次 地 貌 测量 所 得 结果 ,对 一 方形 区 域 (z,y 方向 均 为 从 1 一 10) ,选择 某 些 
地 点 测量 其 相对 于 某 水 平面 高 度 的 数据 ,要 求 用 这 些 数据 尽量 准确 地 绘制 出 该 地 区 的 地 形 
地 狐 。 





表 1-6 某 地 区 地 艇 测 量 结果 













































































1 2 3 生 5 6 7 8 9 10 
二 0.02 一 作 12: 一 和 09 一 0.58 | 一 0.08 
2 0.02 一 2.38 = 一 必 工 
3 0. 10 1.00 一 3.04 一 0.53 0. 10 
于 3.52 
5 一 0.43 | 一 1.98 0.77 2.17 
6 一 2.29 0.69 2.59 0. 30 
莹 一 0.09 | 一 0.31 4.27 一 0.01 
8 5.13 7.40 1. 89 0.04 
9 0.1 0.58 1.75 一 dt 了 
10 一 0.01 0.3 0.01 
解决 此 问题 的 关键 是 要 将 未 测量 地 点 的 高 度 用 数据 搬 值 的 方法 求 出 来 ,然后 就 可 以 用 
MATLAB 绘图 方法 画 出 地 貌 图 。 
编写 以 下 命令 文件 , 即 可 解决 此 问题 。 
x, 中 =meshgrid(1:10); s#% 构 造 测 最 网 格 


h= [0,0.02, -0.12,0,-2.09,0,- 
0.02,0,0, -2.38,0,-4.96,0,0,0 






10,0; 





0,0,- 2.29,0,0.69,0,2.59,0,0.3 


0.58,- 0.08,0,03 


0.13.- 
0,0.1,1,0，- 3.04,0,，- 0.53,0,0.1,05... 


， 工 98,04040,0.7 0,2.17,0,0 天 





-0.09,-0.31,0,0,0,4.27,0,0,07 一 0 
0,0,0%5.13,7.4,0,1.89,0)0.04y035 


下 巧 0,0.58,0 而 , 必 75.0, 二 0.110,0 

0 0.010,00.31000,0,0101， # 测量 数据 点 
[xi,yi=meshgridC1:01;10)) s#% 构造 插值 网 格 
hi = interp2(xyy,hy,xiyyi，spline'); #% 二 维 插值 命令 
Surf(hi)y s% 画 出 地 貌 图 
xlable( x ) ,ylableCY ),zlableC h ); s# 添加 坐标 名 称 


运行 此 文件 ,得 到 如 图 1 - 11 所 示 的 结果 , 即 为 该 地 区 的 地 瑶 形 态 图 。 此 处 用 分 段 三 次 样 


条 插值 会 得 到 更 好 的 效果 。 


1.4.2 夺 灯 光源 投影 区 域 的 绘制 (CUMCM 2002 A) 


2002 年 CUMCM 的 A 题 ,其 中 的 一 个 要 求 是 绘制 投影 区 域 。 需 要 先 建立 车 灯 投 影 的 数 
学 模型 ,然后 再 根据 模型 绘 出 投影 效果 图 。 这 里 不 再 详 述 建 模 过 程 , 仅 给 出 一 段 绘制 投 影 区 域 
的 MATLAB 程序 ,以 说 明 MATLAB 绘图 技术 在 建 模 比 赛 中 的 重要 性 及 程序 绘图 的 方法 。 

利用 MATLAB, 根 据 得 到 的 线 光 源 长 度 , 用 投 点 法 可 画 出 测试 屏 上 反射 光 的 亮 区 ， 


MATLAB 程 序 如 下 : 











图 1-11 地 瑶 图 








p= 0.03;x= 25.0216; 
for L = -0.002:0.0004:0.002 
Y0=(- 0.036:0.001:0.036) * ones(1,73)， 


z0 = ones(73,1) wx(- 0.036:0.001:0.036)4 

x0 = (Y0.-2+z0.-2)/(2x p)5 

5m= (pr3+4xx0xw2xp.wx0+pxw(-4xwyLwyO+3#2xwpxwx0))./(2x(P2+2xwpwx0))3 
yn=(2xpxx0.xwy0+Ep2x#(- 噶 +y0)+YLw(Y0-2-z0.-2))./(pr2+2#pxx0)3 

zn= (pr2+2wpxwx0+2xYL%y0).wz0/(p22wpwx0); 





Y=Y0+(m-yY0).w(x-x0)./(xn-x0)f 
z=z0+(zn-2z0).xw (x-x0)./(xn-x0) 
plot(Y,z，b.) 

xlabelCY ) 

nabel(z ) 

hold on 





即 得 到 反射 光 在 测试 屏 上 的 亮 区 ,如 图 1- 12 所 示 。 
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图 1- 12 反射 光 在 测试 屏 上 的 亮 区 


第 2 章 
| 规划 问题 的 MATLAB 求解 


规划 类 问题 是 常见 的 数学 建 模 问 题 , 离 散 系统 的 优化 问题 一 般 都 可 以 通过 规划 模型 来 求 
解 。 所 以 ,在 建 模 竞赛 中 ,能 够 快速 求解 规划 类 问题 是 数学 建 模 队 员 的 基本 素质 。 利 用 
MATLAB 提 供 的 强大 的 规划 模型 求解 命令 ,可 以 简单 快速 地 得 到 所 要 的 结果 。 一 般 的 标准 
规划 模型 都 可 以 用 这 些 命令 直接 进行 求解 ,本 章 主要 介绍 常见 规划 模型 的 MATLAB 求 解 , 包 
括 线性 规划 、 整 数 规划 和 非 线 性 规划 三 个 部 分 。 掌 握 这 几 部 分 的 操作 ,可 以 解决 大 部 分 的 规划 
模型 的 求解 问题 。 


2.1 线性 规划 








人 们 在 生产 实践 中 ,经 常会 遇 到 如 何 利用 现 有 资源 安排 生 疼 , 以 取得 最 大 经 济 效益 的 问 
题 。 此 类 间 题 构成 了 运筹 学 的 一 个 重要 分 支 一 一 数学 规划 * 而 线性 规划 (linear programming， 
LP) 则 是 数学 规划 的 三 个 重要 分 支 。 历 史上 线性 规划 理论 发 展 的 重要 进程 如 下 

1947 年 美国 数学 家 G、B，Dantzig 提出 线性 规划 的 一 般 数学 模型 和 求解 线性 规划 问题 的 
通用 方法 一 -单纯 形 法 ,为 这 门 学 科 葛 定 了 基础 ; 

1947 年 美国 数学 家 ]. von 诺 伊 曼 提出 对 循 理 论 , 开 创 了 线性 规划 的 许多 新 的 研究 领域 ， 
扩大 了 它 的 应 用 范围 和 解 题 能 力 ; 

1951 年 美国 经 济 学 家 T. C. 库 普 晕 斯 把 线性 规划 应 用 到 经 济 领 域 , 为 此 与 康 托 罗维奇 一 
起 获 1975 年 诺 贝 尔 经 济 学 奖 。 

目前 ,线性 规划 理论 趋向 完善 ,应 用 范围 不 断 延伸 ,已 经 渗透 到 众多 领域 ,特别 是 在 计算 机 
能 处 理 成 千 上 万 个 约束 条 件 和 决策 变量 的 线性 规划 问题 之 后 ,线性 规划 的 适用 领域 更 为 广泛 
了 ,已 成 为 现代 管理 中 经 常 采用 的 基本 方法 之 一 。 


2.1.1 线性 规划 的 实例 与 定义 


【 例 2- 1 中 央 电视 台 为 改版 后 的 非常 6 十 1) 栏 目 播放 两 套 宣传 片 。 其 中 ,宣传 片 甲 播 
肌 时 间 为 3. 5min, 广 告 时 间 为 30s, 收 视 观众 为 60 万 ;宣传 片 乙 播映 时 间 为 lmin, 广 告 时 间 为 
1min, 收 视 观 众 为 20 万 。 广 告 公 司 规定 每 周至 少 有 3. 5 min 广告 ,而 电视 台 每 周 只 能 为 该 栏 
目 宣传 片 提供 不 多 于 16 min 的 节目 时 间 。 电 视 台 每 周 应 播映 两 套 宣传 片 各 多 少 次 ,才能 使 得 
收视 观众 最 多 ? 

分 析 : 建 模 是 解决 线性 规划 问题 的 极为 重要 的 环节 与 技术 。 一 个 正确 数学 模型 的 建立 要 
求 建 模 者 熟悉 规划 问题 的 生产 和 管理 内 容 , 明 确 目 标 要 求 和 错综复杂 的 约束 条 件 。 本 题 首 先 
将 已 知 数据 列 成 表 2 - 1 所 示 的 清单 。 





表 2-1 题 意 信息 清单 














类 别 播放 片 甲 播放 片 乙 节目 要 求 
片 集 时 间 /min 3.5 1 
去 16 
广告 时 间 /min 0.5 1 三 3.5 
收视 观众 /万 人 60 20 























设 电 视 台 每 周 应 播映 片 甲 工 次 , 片 乙 > 次 ,总 收视 观众 为 = 万 人 


目标 函数 ， maxz 一 60z 十 20y (2-1) 
4 十 2y 生 16 

约 东 条件， 5z 十 y 过 3.5 (2- 人 2 
r,yEN 


其 中 ,变量 r,y 称 为 决策 变量 , 式 (2 - 1) 被 称 为 问题 的 目标 函数 , 式 (2 - 2) 中 的 几 个 不 等 
式 是 问题 的 约 东 条 件 , 记 为 st ( 即 subject to)。 上 述 即 为 一 规划 问题 数学 模型 的 三 个 要 素 。 
由 于 上 面 的 目标 函数 及 约束 条 件 均 为 线性 函数 , 故 被 称 为 线性 规划 问题 。 

在 解决 实际 问题 时 ,把 问题 归结 成 一 个 线性 规划 数学 模型 是 很 重要 的 一 步 , 但 往往 也 是 困 
难 的 一 步 ,模型 建立 的 是 否 恰当 直接 影响 到 求解 。 而 选取 适当 的 决策 变量 ,是 建立 有 效 模型 的 
关键 之 一 4 


2.1.2 线性 规划 的 MATLAB 标准 形式 


线性 规划 的 目标 函数 可 以 是 求 最 大 值 ,也 可 以 是 求 最 小 值 :约束 条 件 的 不 等 号 可 以 是 小 于 
号 ,也 可 以 是 大 于 号 。 为 了 避免 形式 多 样 性 带 来 的 不 便 ,MATLAB 中 规定 线性 规划 的 标准 形 
式 为 





miineIzfz st hx 入 1 


其 中 \c 和 x 为 习 维 列 向 量 ;大 为 友 维 列 向 量 ;4 为 四 X7 矩阵 。 
例如 线性 规划 ， 


inaxcTRS. tx 二 六 


的 MATLAB 标准 型 为 


min 一 cx st 一 4xr 安 一 
2.1.3 线性 规划 问题 解 的 概念 


一 般 线 性 规划 问题 的 标准 型 为 
二 (2-3) 





入) 站。 iv25oy 于 (2-4) 


满足 约束 条 件 式 (2- 4) 的 解 xz-Cmn :av) 称 为 线性 规划 问题 的 可 行 解 ,而 使 目标 函数 
式 (2- 3) 达 到 最 小 值 的 可 行 解 称 为 最 优 解 。 所 有 可 行 解构 成 的 集合 称 为 问题 的 可 行 域 , 记 为 尺 。 


2.1.4 求解 线性 规划 的 MATLAB 解法 
自从 G、B. Dantzig 于 1947 年 提出 单纯 形 法 ,60 多 年 来 , 虽 有 许多 变形 体 已 被 开发 ,但 却 


保持 着 同样 的 基本 观念 。 由 于 有 如 下 结论 : 若 线性 规划 问题 有 有 限 最 优 解 , 则 一 定 有 某 个 最 优 
解 是 可 行 区 域 的 一 个 极点 。 基 于 此 ,单纯 形 法 的 基本 思路 是 : 先 找 出 可 行 域 的 一 个 极点 , 据 一 
定 规则 判断 其 是 否 最 优 , 和 否则 转换 到 与 之 相 邻 的 另 一 极点 ,并 使 目标 函数 值 更 优 ; 如 此 下 去 , 直 
到 找到 某 一 最 优 解 为 止 。 这 里 不 再 详细 介绍 单纯 形 法 ,有 兴趣 的 读者 可 以 参看 其 他 有 关 线 性 
规划 的 书籍 。 下 面 介 绍 线性 规划 的 MATLAB 解法 。 
MATLAB 中 线性 规划 的 标准 型 为 
mincTx st 4r 雯 ) 
基本 函数 形式 为 linprog(c,A,b) , 它 的 返回 值 是 向 量 * 的 值 。 还 有 其 他 一 些 函 数 调用 形 
式 (在 MATLAB 指令 窗口 运行 help linprog 可 以 看 到 所 有 的 函数 调用 形式 ) ,如 
[x,fvaI]=linprog(c,A,b,Aeq,beq,LB,UB,X0,OPTIONS) 
其 中 ,fval 返回 目标 函数 的 值 ;Aeq 和 beq 对 应 等 式 约束 hx 一 p;LB 和 UB 分 别 是 变量 x 
的 下 界 和 上 界 ;X0 是 x 的 初始 值 ;OPTIONS 是 控制 参数 。 
【 例 2- 2 求解 线性 规划 问题 : 
min = 一 2z: 十 3zs 十 zs 
ri 十 4rs 十 2r 过 8 
St. 可 十 2zs 之 6 
rz 艺 0 


解 :编写 MATLAB 程序 如 下 : 


c=[2;3;1]; 
a= [1,4,2;3,2,0]; 
b= [8;6]; 
[x, 归 = linprog(c, -a, - b,[],[],zeros(3:D) 
【 例 2-3】〗 求 下 面 的 优化 问题 ; 
imin 司 一 一 5zl 一 4zs 三 6zs 
屋 , Ex: 计 怀 二 罗 
zt 257: 4 总 科 亿 
SN 
3zi 十 2z 委 30 
0 和 mm0 委 ra0 委 mm 
解 : 
DA 
>> 有 


[-55 一 4 -6]; 

下 

>> b = [20; 42; 30]; 

>> lb = zeros(3,1); 

>> [xfval,exitflagv,output,lambda] = linprog(f,ab'[],[],1b) 


结果 为 


[0 


x = s# 最 优 解 
0.0000 
15.0000 
3.0000 

fwal = s# 最 优 值 





一 78.0000 


exitflag = 多 收 敛 
二 
output = 
iterations: 6 当选 代 次 数 
cgiterations: 0 
algorithm: lipsol 多 所 使 用 规则 
lanbda = 
ineqlin: [3xl double] 





eqlin; [oxl double] 
upper: [3xl double] 
lower: [3xl double] 
>> lanmbda. ineqlin 
ans = 
0.0000 
1.5000 
0.5000 
>> lambda, lower 
ans = 
1.0000 
0.0000 
0.0000 


【 例 2=4】 求解 下 列 线性 规划 问题 : 
maxz 三 2zr 十 37 二 5zs 
7 十 mHz 一 7 
st 12 522 十 7 志 10 
Ravzrs 三 0 
解 ; 
(1) 编写 M 文件 。 


Re [1385 

a=[-2,5,-1];b= -10; 

aeq= [1,1,1]; 

beq=7; 

% 是 求 最 大 值 而 不 是 最 小 值 , 注 意 这 里 是 "- c" 而 不 是 "c” 
Xx= linprog( - cyaybyaeqybeq,zeros(3,1)) 

value= cx 


(2) 将 M 文件 存盘 ,并 命名 为 examplel. m。 

(3) 在 MATLAB 指令 窗口 运行 examplel 即 可 得 所 求 结 果 。 

【 例 2 - 5】 求解 下 列 线性 规划 问题 

一 170. 858 2z, 一 17.725 4zs 十 41.258 2zrs 十 2. 218 2zr4 十 131. 818 2zs 一 500 000 





maxzx 





一 0.17037zr, 一 0.5324zs 十 zs 妇 0 
0.17037zs 十 0.53247s 二 888 115 
痉 十 32%zs 十 zs 委 166805 
st ra 委 521265. 625 
2 十 zi 过 683400 
2 十 zs 660000 
和 0 二 8 二 
解 :为 了 便于 求解 ,将 上 述 求解 最 大 值 线性 规划 问题 转化 成 求解 最 小 值 问题 : 


minx“ 一 一 170. 858 2z: 十 17.725 47s 一 41. 258 2zs 一 2. 218 2z, 一 131. 818 2zs 十 500 0007re 


:一 0.17037zrs* 一 0.5324zs 十 zs 委 0 
0.17037zs 十 0.5324z; 二 888115 
十 32%zs 十 zs 过 166805 
zs 二 521265. 625 
2 十 Ti s683400 
一 一 zs 魏 一 660000 
去 二 1 

0 一 站 3 





MATLAB 源 程序 : 


上 =[-170.8582 17.7254 - 41.2582 - 2.2182 -131.8182 500000]; 
R= [1 -0.17037 -0.5324010;0 0.17037 0.5324000;10.321000;010000;001100;000 


二 2 


即 


b= [0;888115;166805;521265. 625;683400; -= 660000]; 
aeq=[000001] 

beq= [1]; 

lb = 下 0 不 0:0:00， 内 

[xifval 2 - 1 ipEog(f 银 汪 ,Meq,beg, 卫 ,[]) 


程序 输出 结果 


区 = 
1.0e+ 005 关 
:00000000000000 
:70617739889132 
.12207323235472 
.71192676764526 
.88807323235476 
. 00001000000000 
fwal = -1.407864558820066e+ 007 


0 
I 
1 
上 
0 
0 


ZI 一 0， 272 一 170618， zs 二 112207， 一 571193， zs 一 88 807， ze 一 1 
minz' 一 一 14078 646， maxz 一 一 minz' 一 一 (一 14078646) 一 14078 646 





2.2 斐 线 性 规划 


2.2.1 非 线性 规划 的 实例 与 定义 


如 果 目 标 函 数 或 约束 条 件 中 包含 非 线性 函数 ,就 称 这 种 规划 问题 为 非 线性 规划 问题 。 一 
般 说 来 , 解 非 线性 规划 要 比 解 线性 规划 问题 困难 得 多 ,而 且 非 线性 规划 也 不 像 线性 规划 那样 有 
单纯 形 法 这 一 通用 方法 。 非 线性 规划 目前 还 没有 适 于 各 种 问题 的 一 般 算法 ,各 种 方法 都 有 自 
己 特 定 的 适用 范围。 
下 面 通过 实例 归纳 出 非 线性 规划 数学 模型 的 一 般 形式 ,介绍 有 关 非 线性 规划 的 基本 概念 。 
【 例 2 - 6】 投资 决策 问题 : 某 企 业 有 ， 个 项 目 可 供 选择 投资 ,并 且 至 少 要 对 其 中 一 个 项 
目 投资 。 已 知 该 企业 拥有 总 资金 A 元 ,投资 于 第 1G 一 1,2,…,z) 个 项 目 需 花 资金 w 元 ,并 杞 
计 可 收益 元 。 试 选择 最 佳 投资 方案 。 
解 : 设 投资 决策 变量 为 ; 
，_ | 1， 决定 投资 第 ; 个 项 目 
” | 0， 决定 不 投资 第 去 个 项 目 ” 


则 投资 总 额 为 ai 投资 总 收益 为 > 5 。 
因为 该 公司 至 少 要 对 一个 项 目 投资 ,并 且 总 的 投资 金额 不 能 超过 总 资 使 责 , 故 有 限制 条 件 ， 
0 员 


由 于 zi Gi 一 1,2,…,) 只 取 值 0 或 1 记 以 还 有 
了 QI 二 0，1 一 1,2，on 
另外 ,该 公司 至 少 要 对 一 个 项 绅 投 资 , 因 此 有 
澡 卡 本 
最 佳 投资 方案 应 是 投资 额 最 小 而 总 收益 最 大 的 方案 ,所 以 这 个 最 佳 投资 决策 问题 归结 为 


在 总 资金 以 及 决策 变量 ( 取 0 或 1) 的 限制 条 件 下 , 极 大 化 利润 ( 即 总 收益 ) 与 总 投资 之 差 。 因 


1 2 


maxQ 一 下 共产 光 Ti 
= :1 


补 =>1 
例 2- 6 是 在 一 组 等 式 或 不 等 式 的 约束 下 , 求 一 个 函数 的 最 大 值 (或 最 小 值 ) 问 题 ,其 中 目 
标 函 数 或 约束 条 件 中 至 少 有 一 个 非 线性 函数 ,这 类 问题 称 之 为 非 线性 规划 问题 , 简 记 为 (CNP) 。 
可 概括 为 一 般 形 式 : 
minfCxr) 





电 七 而) 过 0 天 一 1 2v…d 
机 人 一 而 二 一生 ee 访 
其 中 ,zx 一 [ri xz … 本 称 为 模型 (NP) 的 决策 变量 ;j 称 为 目标 函数 ;g&: G 一 1， 
2,…，, 记 和 太 G 一 1,2,…,q) 称 为 约束 函数 。 另 外 ,g;(Cz) 一 0G 一 1,2,…,/) 称 为 等 式 约 东 ; 
态 Cr) 过 0G 王 1,2,…,9) 称 为 不 等 式 约束 。 


2.2.2 非 线性 规划 的 MATLAB 解法 
在 MATLAB 中 非 线性 规划 的 数学 模型 可 写成 以 下 形式 ， 


min Fr) 
4xr 委 了 8 
Aeq.x 王 Beq 
“|]co)y 和 0 
CeqGCr) 一 0 
其 中 ,jz) 是 标量 函数 ;4 ,B,Aeq,Beq 是 相应 维 数 的 矩阵 和 向 量 ;CCx),Ceq(x) 是 非 线性 
向 量 函 数 。 
MATLAB 中 的 命令 是 


X = FMINCONCFUN,X0 ,A,B,Req,Beq,LB,UB,NONLCON,OPTIONS) 


它 的 返回 值 是 向 量 x。 其 中 ,FUN 是 用 M 文 件 定义 的 函数 fCz) 。X0 是 玖 的 初始 值 。A,B， 
Aeq,Beq 定 义 了 线性 约束 hx 扩 B,Aeq . x 王 Beq, 如 果 没 有 线性 约束 , 则 入 ==[],B=[,Aeq 一 
[ 口 ,Beq=[。LB 和 TUB 是 变量 x 的 下 界 和 上 界 * 如 果 上 界 和 下 界 没有 约束 , 则 LB=[],UB= 
[ 品 ; 如 果 z 无 下 界 , 则 LB= 一 inf; 如 果 * 无 上 界 , 则 UB=inf。NONLCON 是 用 M 文件 定义 
的 非 线 性 向 量 函 数 CC(x) ,CeqCx)。OPTIONS 定义 了 优化 参数 ,可 以 使 用 MATLAB 默认 的 
参数 设置 。 
【 例 2=7] 求 下 列 非 线 性 规划 问题 ; 
minrGa) 三 巡 十 于 十 8 


一 za 人 0 
志 4 一 站 一 二 2 二 9 


zsza 30 
(1) 编 写 M 文 件 funl. m 和 M 文件 fun2. m。 


function f = funl(x) 
上 = x(1)“2+x(2)-2+83 


function [g,bh] = fun2(x)5 
g= 一 x(1)"2+x(2)5 
h= -x(1) - x(2)-2+2; #% 等 式 约束 


(2) 在 MATLAB 的 命令 窗口 依次 输入 如 下 语句 : 
options = optimset 


[x，, = 甸 incon( funl ,rand(2,1),[, 口 , 口 , 口 ,zeros(2,1) 口 ，..- 
“fun2 ，options) 


就 可 以 求 得 当 六 一 1,z: 一 1 时 ,最 小 值 * 一 10。 

【 例 2- 8】 求 下 列 非 线性 规划 问题 : 
maxz= 一 Vzi 十 Vzs 十 Vzs 十 Vzte 
zi <400 
1.1zr 十 zs 和 440 

s.t. 11.21z: 十 1.1zrs 十 zs 近 484 

1.331zi: 十 1.21zs 十 1.1zs 十 zi 委 532.4 
mi 全 0，i 一 1,2,3,4 

01) 编写 M 文件 ,定义 目标 函数 : 





function 下 = fun44(x) 
于 = - (sqrt(x(1)) + sqrt(x(2)) + sqrt(x(3)) + sqrt(x(4))); 


(2) 编 写 M 文件 ,定义 约束 条 件 : 


function [g,ceq] = myconl(x) 

g(1) =x(1) -400; 

.1wx(1) +x(2) -440; 

.21wx(1)+1.1wx(2)+x(3) -4844 

.331*wx(1) +1.21wx(2) +1.1%x(3) +x(4) -532.43 





《3) 编写 主 程序 ( 既 可 以 编写 M 文件 ,也 可 以 在 Command Window 直接 输入 命令 ) , 


x0= [1;1;1;1];lb= [0;0;0;0]iub=[]iR= 口 jp=[];Meq=[];beq=[]; 
[x,fval] = fmincon( fun44 ,x0,R,b,heqybeq,lbyub， myconl ) 


程序 输出 结果 为 


x = 
1.0e+ 002 * 
0. 84243824470856 
1.07635203745600 
1.28903186524063 
1.48239367919807 

fwal = 

一 43.08209516098581 


所 以 最 终结 果 为 
[ 一 84.24， za 一 107.63， 2 一 128.90， zt 一 148.23 
| = = 一 (一 43.08) 一 43.08 


2.2.3 二 次 规划 

若 某 非 线性 规划 的 目标 函数 为 自 变量 z 的 二 次 函数 ,约束 条 件 又 全 是 线性 的 , 则 称 这 种 
规划 为 二 次 规划 。 

在 MATLAB 中 二 次 规划 的 数学 模型 可 表述 如 下 : 





min 到 二 rz，st hr 二 


其 中 ,1 是 列 向 量 ;4 是 相应 维 数 的 矩阵 ; 瑟 是 实 对 称 矩 阵 。 实 对 称 和 矩阵 的 定义 :如 果 有 
导 阶 矩阵 4 ,其 各 个 元 素 都 是 实数 , 且 满足 cy 一 au ( 转 置 为 其 本 身 ), 则 称 4 为 实 对 称 和 矩阵。 
MATLAB 中 求解 二 次 规划 的 命令 是 


[X,FVAL] = QUADPROG(H, 上 ,Rb,heqybeq,LB,UB,X0,OPTIONS) 


X 的 返回 值 是 向 量 *x;FVAL 的 返回 值 是 目标 函数 在 * 处 的 值 。 具 体 细 节 可 以 参看 在 
MATLAB 指令 中 运行 help quadprog 后 的 帮助 。 
【 例 2-9】 求解 二 次 规划 : 
minfx) 一 2 寻 一 4zuzs 十 4 太一 6zl 一 3zz 


zi 十 zs 和 3 
st 14z: 十 zz 委 9 
lszs 之 0 
解 :编写 如 下 程序 : 
h=[4,-4:-4,8]; 
f=[-6-3]? 
a=[1,1;4,1]; 
b= [3;9]， 


[x,value] = quadprog(h,f,a,b,[],[,zeros(2,1)) 


求 得 
1.9500T 
7 一 位 minjfCxr) 一 一 11.0250 

利用 罚 丽 数 法 ,可 将 非 线 性 规划 问题 的 求解 转化 为 求解 一 系列 无 约束 极 值 问题 ,因而 也 称 
这 种 方法 为 序列 无 约束 最 小 化 技术 , 简 记 为 .SUMTI(sequential unconstfained minimization 
technique) 。 

罚 函 数 法 求解 非 线 性 规划 问题 的 思想 是 ,利用 问题 中 的 约束 函数 做 出 适当 的 罚 函 数 ,由 此 
构造 出 带 参数 的 增 广 目标 函数 ,把 问题 转化 为 无 约束 非 线性 规划 问题 。 罚 函数 主要 有 两 种 形 
式 , 一 种 称 外 罚 本 数 法 ; 另 一 种 称 内 罚 杖 数 法 。 下 面 主要 介绍 外 罚 函 数 法 。 

考虑 如 下 问题 ; 

minf Co) 
网 检 二 病人 
8 二 人 ii(z) 世 0， 宇 一 12s 
和 Air) 一 0， 守 一 1,2.…… 
取 一 个 充分 大 的 数 M 二 0, 构 造 函 数 


PCr,AMD) = fr) 十 M>)max(gi(r),0) 一 M>)min(u(Cz),0) 二 My | 大 Cz) | 
当 气 评 


或 PCxr, MD 一 帮 (xz) 十 Mimax(GCr) ,0) 十 MaminCHCxz),0) 十 Ms | KCxz) | 
其 中 ， 


Cr) Cz) 太 Cz) 


号 和 而 本 
训 寺 | 本 全 | 大 一 | 全 有 |， 训 亲 二 后 
机 | 太 CD) 


Mi ,Ma ,Ms 为 适当 的 行 向 量 ,MATLAB 中 可 以 直接 利用 max 和 min 函数 。 则 以 增 广 目标 郴 
数 PCx,M) 为 目标 函数 的 无 约束 极 值 问题 
minPCr,MD) 
的 最 优 解 x 也 是 原 问题 的 最 优 解 。 
【 例 2- 10】 求 下 列 非 线性 规划 : 
minf xz) 一 好 十 肆 十 8 


一 zz 过 0 
st 1 一 一 妇 十 2 一 0 


riyza 之 0. 
解 : 
(1) 编写 M 文件 test m。 


function g = test(x) 

M= 500001 

于 = x(1)"2+x(2)-2+85 

g=f-Mxmin(x(1),0) -Mw min(x(2),0) -NMwminCx(I)-2 一 x(2),0).，. 
+M% abs(-x(1) 一 x(2)"2+2)3 


(2) 在 MATLAB 命令 窗口 输入 : 
[x, 轨 = fminunc( test' ,rand(2,1)) 


即 可 求 得 问题 的 解 。 


2.3 整改 规划 


2.3.1 整数 规划 的 定义 


规划 中 的 变量 (部 分 或 全 部 ) 限 制 为 整数 时 , 称 为 整数 规划 。 若 在 线性 规划 模型 中 ,变量 限 
制 为 整数 , 则 称 为 整数 线性 规划 。 目 前 所 流行 的 求解 整数 规划 的 方法 ,往往 只 适用 于 整数 线性 
规划 。 目 前 还 没有 一 种 方法 能 有 效 地 求解 一 切 整数 规划 。 常 见 的 整数 规划 问题 的 求解 算法 有 
以 下 几 种 ， 

(1) 分 枝 定 界 法 :可 求 纯 或 混合 整数 线性 规划 。 

(2) 制 平 面 法 :可 求 纯 或 混合 整数 线性 规划 。 

53) 隐 枚 举 法 :用 于 求解 0- 1 整数 规划 ,有 过 滤 隐 枚 举 法 和 分 枝 隐 枚 举 法 。 

(4) 匈牙利 法 :解决 指派 问题 (0 - 1 规划 特殊 情形 ) 。 

(5) 蒙特 卡 罗 法 :求解 各 种 类 型 规划 。 


2.3.2 0-1 整 数 规划 


0 -1 整数 规划 是 整数 规划 中 的 特殊 情形 , 它 的 变量 *; 仅 取 值 0 或 1, 这 时 z 称 为 0-1 变 
量 ,或 称 二 进 制 变量 。z 仅 取 值 0 或 1, 这 个 条 件 可 由 下 述 约束 条 件 表示 :0 入 和 1,m EN 或 
2 一 2) 一 0 一 1,2,… 和 2。 在 实际 问题 中 ,如 果 引入 0-1 变量 ,就 可 以 把 有 各 种 情况 需要 分 
别 讨 论 的 线性 规划 问题 统一 在 一 个 问题 中 讨论 了 。 
下 面 举 例 说 明 一 种 解 0- 1 整数 规划 的 隐 枚 举 法 。 所 谓 隐 枚 举 法 是 指 在 枚 举 过 程 中 引信 
过 滤 条 件 ,而 与 之 对 应 的 是 显 枚 举 法 (一 种 穷 举 法 ) 。 
【 例 2- 11】 max 一 3z1 一 2rs 十 5zs 
十 2m 一 委 2 
Za 十 4z 十 ms 委 4 
s.t. 1z 十 zz 委 3 
4r 十 zs 和 6 
zivzanzs 一 0 或 1 
求解 思路 及 改进 措施 : 
(1) 先 试探 性 地 求 一 个 可 行 解 , 易 看 出 (ziyza52zs)= (170307 满 足 约束 条 件 , 故 为 一 个 可 
行 解 , 且 相 应 的 目标 函数 值 为 = 一 3。 
(2) 因为 是 求 极 大 值 问题 , 故 求 最 优 解 时 ,凡是 目标 值 >< 王 3 的 解 不 必 检验 是 否 满足 约束 
条 件 即 可 删除 ; 因 它 肯定 不 是 最 优 解 , 于 是 应 增加 一 个 约束 条 件 ( 目 标 值 下 界 ) 
3 一 222 十 5z, 亏 3 
称 该 条 件 为 过 滤 条 件 (filtering contraint)。 从 而 原 问题 等 价 于 


max 之 三 -371 一 22 十 573 


3xi 一 2rz 十 5z: 3 (人 -5 
| kz -5b) 
iz 目 十 4 了 目 zy <S4 (中 5@ 
We ( 强 5 
本 (2 -5e) 
5 二 二 磺 浊 C2-5D 


若 用 全 部 枚 举 法 ,3 个 变量 共有 8 种 可 能 的 组 合 ,依次 检验 这 8 种 组 合 , 看 其 是 否 满足 条 
件 式 (2- 5a) 一 式 (2 - 5e), 对 某 个 组 合 , 若 它 不 满足 式 (2 - 5a), 即 不 满足 过 滤 条 件 , 则 
式 (2- 5b) 一 式 (2 - 5e) 即 可 行 性 条 件 不 必 再 检验 ; 若 它 满足 式 (2 - 5a) 一 式 (2 - 5e) 且 相应 的 
目标 值 严 格 大 于 3, 则 进行 步骤 (3) 。 

(3) 改进 过 滤 条 件 。 

(4) 由 于 对 每 个 组 合 首先 计算 目标 值 以 验证 过 滤 条 件 , 故 应 优先 计算 目标 值 >, 这 样 可 提 
前 抬 高 过 滤 门 槛 ,以 减少 计算 量 。 

按 上 述 思路 与 方法 , 例 2- 11 的 求解 过 程 可 由 表 2- 2 表示 。 





表 2-2 例 2-11 的 求解 过 程 






























































(rayzrzyza) | 目标 值 省 全 和 过 旭 条 件 
式 (2-5a) | 式 (2-5b) | 式 (2-5c) | 式 (2-5d) | 式 (2- 5e) 

(0,0,0) 0 X 

《1,0,0) 3 V 忆 V V V 3zl 一 2rz 十 53 关 3 

(0,1,0) 一 2 X 

(0,0,1) 5 V 4 V 以 V 3zl 一 2rz 十 53 六 5 

(1,1,0) 1 X 

《10,1) 8 V V V V V 3zl 一 2rz 十 5 六 8 

(1,1,1) 6 X 

CO,1D 3 尝 

从 而 得 最 优 解 (zy ,zz ,zy ) 一 (1,0,1) ,最 优 值 >" 一 8。 


2.3.3 随机 取样 计算 法 


前 面 介绍 的 常用 的 整数 规划 求解 方法 ,主要 是 针对 线性 整数 规划 而 言 ,而 对 于 非 线性 整数 
规划 目前 尚未 有 一 种 成 熟 而 有 效 的 求解 方法 ,因为 非 线性 规划 本 身 的 通用 有 效 解法 尚未 找到 ， 
更 何况 是 非 线性 整数 规划 。 

然而 ,尽管 整数 规划 由 于 限制 变量 为 整数 而 增加 了 难度 ,然而 又 由 于 整数 解 是 有 限 个 ,于 
是 为 枚 举 法 提供 了 方便 .当然 ;在 自 变量 维 数 很 大 和 取 值 范围 很 宽 的 情况 下 * 企 图 用 显 枚 举 法 
( 即 穷 举 法 ) 计 算出 最 优 值 是 不 现实 的 ,但 是 应 用 概率 理论 可 以 证 明 ,在 有 一 定 计算 量 的 情况 
下 ,完全 可 以 得 出 一 个 满意 解 。 


【 例 2- 12】 


TPR 十 证 笠 直 < 400 


st. 17l 十 2zs 十 2zs 十 zi 十 6zs 魏 800 


2zr: 十 z 十 6zs 200 
2 十 zi 十 5rs 委 200 


2ms 


已 知 非 线性 整数 规划 为 
maxz 们 型 + 十 Se 汪 4zg<2z 一 8zog 一 2 3zs Pa 一 
Q@ 和 如 世 99 有 区 一 1 呈 ,5 


解 :(1) 首 先 编写 M 文件 mente. m, 定 义 目 标 函 数 太 和 约束 向 量 函 数 g ,程序 如 下 : 


function [f,g] = mengte(x) 
上 = x(1)-2+x(2)-2+3*x(3)-2+4xx(4)-2+2xx(5) -8Bxx(1) -2xx(2) 一 3xx(3)-x(4) 一 2xx(5)3 


g(1) = sum(x) - 400; 
g(2) =x(1) +2xx(2)+2#*x(3)+x(4)+6 关 X(C5) 一 8003 


g(3) =2*x(1)+x(2)+6#x(3) 一 200; 
g(4) =x(3) +Xx(4) +5*xXx(5) 一 2003 


(2) 编写 如 下 程序 求 问题 的 解 : 


rand( state ,sum(clock))3 


p0=0; 





tic 
for =1:105 
x= 99x rand(5,1); 
Xl = floor(x) ;x2 = ceil(x)5 
[f,g] = mengte(xl); 
证 sum(g <= 0) ==4 
证 p0 <= 芋 
x0= 了 ipP0=f 
end 
end 
[f,g] = mengteCx2) 
话 sum(g <= 0) == 4 


证 p0 <= 开 
x0=x2;p0= 全 
end 
end 
end 
x0，p0 
toc 


运行 程序 , 即 得 到 如 下 的 输出 结果 : 


44 94 中影 禾 了 和 
47789 

需要 说 明 的 是 ,由 于 是 随机 取样 ,所 以 每 次 运行 的 结果 都 可 能 不 一 致 ,但 误差 不 会 太 大 。 

关于 整数 规划 模型 的 求解 ,虽然 MATLAB 可 以 解决 ,但 对 于 比较 标准 的 整数 规划 模型 ,建议 


还 是 首先 考虑 用 Lingo 求 解 ,这 方面 Jingo. 比 较 有 优势 ,实在 不 行 ,再 考虑 用 MATLAB 来 解 
决 ,MATLAB 的 通用 性 比较 强 。 


第 3 章 
灰色 预测 及 其 MATLAB 实现 


在 数学 建 模 中 经 常会 遇 到 数据 的 预测 问题 ,有 些 赛 题 中 ,数据 的 预测 占 主导 地 位 ,如 表 3- 1 


























所 列 。 
表 3-1 历届 CUMCM 数据 预测 题目 
年 度 类 别 题 目 命题 人 
2003 A 题 SARS 的 传播 问题 CUMCM 组 委 会 
2005 人 题 长 江水 质 的 评价 和 预测 问题 昔 中 庚 
2006 B 题 艾滋 病 疗法 的 评价 及 疗效 的 预测 问题 边 答 萍 
2007 A 题 中 国人 口 增 长 预测 问题 唐 云 

















有 些 问题 则 是 需要 在 求解 的 过 程 中 进行 数据 预测 ,如 .2009 年 CUMCM D 题 "会议 筹 备 ” 
对 与 会 人 数 的 确定 等 。 灰 色 模型 gray model, 又 称 灰色 理论 ) 有 严格 的 理论 基础 ,最 大 优点 是 
实用 。 用 灰色 模型 预测 的 结果 比较 稳定 ,不 仅 适用 于 大 数据 量 的 预测 ,在 数据 量 较 消 时 (数据 
只 要 多 于 3 个 即 可 ) 预 测 结果 依然 较 准确 。 


3.1 记 色 预测 基础 知识 


灰色 系统 理论 认为 :系统 的 行为 现象 尽管 是 有 瞳 胱 的 s 数 据 是 复杂 的 ,但 它 毕竟 是 有 序 的 ,是 
有 整体 功能 的 。 在 建立 灰色 预测 模型 之 前 ' 需 先 对 原始 时 间 序 列 进行 数据 处 理 ,经 过 数据 预 处 
理 后 的 数据 序列 称 为 生成 列 。 对 原始 数据 进行 预 处 理 , 不 是 寻求 它 的 统计 规律 和 概率 分 布 ,而 
是 将 杂乱 无 章 的 原始 数据 列 通过 一 定 的 方法 处 理 , 变 成 有 规律 的 时 间 序 列 数据 , 即 以 数 找 数 的 
规律 ,再 建立 动态 模型 。 灰 色 系统 常用 的 数据 处 理 方式 有 累加 和 累 减 两 种 ,通常 用 累加 方法 。 

灰色 预测 通过 鉴别 系统 因素 之 间 发 展 趋势 的 相 异 程度 ,并 对 原始 数据 进行 生成 处 理 来 寻 
找 系统 变动 的 规律 ,生成 有 较 强 规律 性 的 数据 序列 ,然后 建立 相应 的 微分 方程 模型 ,从 而 预测 
事物 的 未 来 发 展 趋势 。 灰 色 预 测 的 数据 是 通过 生成 数据 的 模型 所 得 到 的 预测 值 的 着 处 理 结 
果 。 灰 色 预 测 是 以 灰色 模型 为 基础 的 ,在 诸多 的 灰色 模型 中 ,以 灰色 系统 中 单 序列 一 阶 线性 微 
分 方程 模型 GM(1,1) 模 型 最 为 常用 。 下 面 简要 地 介绍 GM(1,1) 模 型 。 

设 有 原始 数据 列 x 一 Cr (1) ,zz (2)，…z 9 CD))，m 为 数据 个 数 。 

如 果 根 据 *"” 数据 列 建立 GM(1,1) 来 实现 预测 功能 , 则 基本 步骤 如 下 : 

(1) 原始 数据 累加 以 便 弱化 随机 序列 的 波动 性 和 随机 性 ,得 到 新 数据 序列 : 

xz 一 (rz@(1)vzw(2)5az0 (2)) 

其 中 ,z? (中 各 数据 表示 对 应 前 几 项 数据 的 累加 。 





二 
各 





或 xz 十 1D) 一 -"， 上 一 1,2，… 季 


(2) 对 z9 CD) 建 立 ri) 的 一 阶 线性 微分 方程， 


oa 吕 
+arw 一 x 


其 中 ,as 证 和 有 的 有 效 区 间 是 〈 一 2，2) ,并 记 
4 村 成 的 和 阵 为 5 一 人 1 只 要 求 出 参数 avv, 就 能 求 出 xz (:) ,进而 求 出 x'o 的 未 来 预 
测 值 。 

(3) 对 累加 生成 数据 做 均值 生成 吾 与 常数 项 向 量 Y., 即 
0.5Czo (1) 十 zo (2)) 
0.5Czm (2) 十 ze (3)) | 了 一 Cr (2)5z (3) 和 rz (GD 
0.5Czo CO 一 1) 十 zo Ca)) 
〈4) 用 最 小 二 乘法 求解 灰 参 数 5, 则 


有 一 





2 = 世 re CBTB)-1BTY.， 


中 二 这 进 行 求解 ,得 





(5) 将 灰 参 数 2 代 入 一 +azw 一 


人 计 人 


由 于 @ 是 通过 最 小 三 乘法 求 出 的 近似 值 ,所 以 f (十 1) 是 一 个 近似 表达 式 ,为 了 与 原 序列 
zt 十 1) 区 分 开 来 , 故 记 为 二 2 (t 十 1) 。 

(6) 对 函数 表达 式 F9 (十 1) 及 EC 进行 离散 ,并 将 二 者 做 差 以 便 还 原 *" 原 序列 ,得 
到 近似 数据 序列 大 "(t 十 1D) 如 下 : 

TOTUNRD 0 CPU) 王 Eco) 
《7)》 对 建立 的 灰色 模型 进行 检验 ,步骤 如 下 : 
人 计算 与” (0 之 间 的 残 差 ”( 厅 和 相对 误差 gz): 
e 吧 (Ci) 一 工 中 一 下 加 (t) 


et) 
0] 


4Cr) 一 
@@ 求 原始 数据 *'" 的 均值 以 及 方差 。 
图 求 ” (0) 的 平均 值 7 以 及 残 差 的 方差 ': 。 
图 计算 方差 比 C= 汪 。 
图 求 小 误差 概率 P 一 P{|e(z)| 一 0.6745s:} 。 
@@ 灰色 模型 精度 检验 如 表 3- 2 所 列 。 
表 3- 2 灰色 模型 精度 检验 对 照 表 





等 级 相对 误差 4 方差 比 C 小 误差 概率 已 
工 级 去 0.01 天 0.35 >0.95 








下 级 < 二 0.05 到 0.50 < 一 0.80 























续 表 3-2 











等 级 相对 误差 9 方差 比 C 小 误差 概率 也 
卫 级 一 0.10 < 到 0.65 < 到 0.70 
级 盖 0.20 盖 0.80 < 到 0.60 




















在 实际 应 用 过 程 中 ,检验 模型 精度 的 方法 并 不 唯一 。 可 以 利用 上 述 方法 进行 模型 的 检验 ， 
也 可 以 根据 vd(z) 的 误差 百分比 并 结合 预测 数据 与 实际 数据 之 间 的 测试 结果 酌情 认定 模型 是 
否 合理 。 

(8) 利用 模型 进行 预测 : 


生 册 色 《17 5 元史 (27 人 人 (71 天 四 (十 1) 和 全 oa| 





原 数列 的 模拟 未 来 数列 的 预测 


3.2 记 色 预测 的 MATLAB 程序 


3.2.1 呐 型 程序 结构 


灰色 预测 中 有 很 多 关于 矩阵 的 运算 ,这 可 是 MATLAB 的 特长 ,所 以 用 MATLAB 是 实现 
灰色 预测 过 程 的 首选 。 用 MATLAB 编写 灰色 预测 程序 时 ,可 以 完全 按照 预测 模型 的 求解 步 
了 又, 即 

《1) 对 原始 数据 进行 累加 。 

(2) 构造 累加 矩阵 瑟 与 常数 向 量 。 

(3) 求解 灰 参 数 。 

(4) 将 参数 带 和 人 预测 模型 进行 数据 预测 。 

下 面 以 某 公 司 收入 预测 问题 为 例 介绍 灰色 预测 的 MATLAB 实现 过 程 。 

已 知 某 公 司 1999 一 2008 年 的 利润 为 (单位 :元 /年 ):[89677,99 215, 109 655,120 333， 
135823,159 878,182321,209407,246 619,300 670], 现 在 要 预测 该 公司 未 来 几 年 的 利润 情况 。 

具体 的 MATLAB 程序 见 P3 - 1。 























程序 编号 | P3- 1 | 文件 名 称 | main0301. m | 说 明 灰色 预测 公司 的 利润 
clear 
syms a bi 
c=[abJ; 
有 = [89677,99215,109655,120333,135823,159878,182321,209407,246619,300670]; 
B= cumsum(a); # 原 始 数据 累加 


n= lengthCGBD; 
for i=1:Cn-1) 
CCGiD = (BCGi) + BCGE+1))/2; 多 生 成 楷 加 矩阵 
end 
s* 计算 待定 参数 的 值 
D=aiDGD) =[]; 
D=D， 
了 =[- Ciones(1,n-1)]; 














C= inv(Ex 卫 ) 关 瑟 #D5 
c=c3 
a=c(1);b=c(2)5 
# 预测 后 续 数 据 
?= 口 ;CD = MI)3 
for i=2:Cn+10) 
FGiD = (MI) - b/a)/exp(a* (i 1))+b/ay 
end 
SG=[];GGD) = MD 
for i=2:Cn+10) 
GGi) =FCGi) -了 GE 一 1)3 s# 得 到 预测 出 来 的 数据 
end 
tl = 1999:2008; 
t2 = 1999:2018 
G 
plot(t1,R,7o",t2,G) s# 原始 数据 与 预测 数据 的 比较 











运行 该 程序 ,得 到 的 预测 数据 如 下 : 


G = 
1.0e+ 006 * 
Columns 1 through 14 
0.0897 ”0.0893 ”0.1034 ”0.1196 ”0.1385 0.1602 
0.2146 0.2483 0.2873 0.3325 0.3847 、 0.4452 0.5152 
Columns 15 through 20 
0.5962 0.6899 ”0.7984 。 0.9239 “1:-0691 。 1.2371 


该 程序 还 显示 了 预测 数据 与 原始 数据 的 比较 图 ,如 图 3- 1 所 示 。 


四 





利 臣 下 





陋 5 5 
皇 份 放 


图 3-1 某 公司 利 润 预测 数据 与 原始 数据 的 比较 


3.2.2 灰色 预测 程序 说 明 
(D 先 熟 悉 程序 中 各 条 命令 的 功能 ,以 加 深 对 灰色 预测 理论 的 理解 。 


0.1854 


(2) 在 实际 使 用 时 ,可 以 直接 套用 该 段 程序 ,把 原 数据 和 时 间 序 列 数据 替换 就 可 以 了 。 





者 也 可 以 对 预测 数据 进行 其 他 方式 的 精度 检验 。 


3.3 友 色 预测 应 用 实例 


(3) 模型 的 误差 检验 可 以 灵活 处 理 ,图 3- !1 给 出 的 是 预测 数据 与 原始 数据 的 比较 图 , 读 








3.3.1 实例 一 :长 江水 质 的 预测 (CUMCM 2005 A) 

长 江 的 水 质问 题 是 一 个 复杂 的 非 线性 系统 ,但 是 由 于 数据 样本 少 ,需要 预测 的 时 间 长 , 直 
接应 用 神经 网 络 很 难 取得 理想 的 效果 。 考 虑 到 污水 排放 量 的 变化 规律 是 一 个 不 确定 的 系统 ， 
且 本 题 给 出 污水 排放 量 数据 样本 比较 少 ,还 要 求 做 出 长 达 10 年 的 预测 ,因此 采用 灰色 预测 方 
法 来 预测 未 来 的 污水 排放 量 。 

对 原 题 附件 4 中 的 数据 进行 整理 可 以 得 到 10 年 的 长 江 污水 量 排放 数据 ,如 表 3- 3 所 列 。 

表 3-3 1995 一 2004 年 长 江 污水 排放 量 
年 份 1995 | 1996 | 1997 | 1998 | 1999 | 2000 | 20ol 2002 | 2003 | 2004 
污水 基 / 亿 吨 | 174 179 183 189 207 234_ |.220.5 卜 256 270 285 



































以 P3 = 1 的 程序 段 为 基础 ,将 表 3 - 3 的 数据 带 和 ,并 更 新 时 间 轴 数据 ， 


程序 输出 的 图 形 如 图 3=- 2 所 示 。 








即 得 到 新 程序 P3- 2。 





程序 编号 | P3 - 2 | 文件 名 称 





Changjiang.mm 


说 明 





预计 长 江 污水 





放量 





clear 
syma a bi 
c=[abJ， 
RM=[174 179 
B= cumsun(CR); 
nz= lengthCh); 
fori=i:tn-1) 
Ci) =(B(i) 十 B( 主 +177/25 


183 189 


end 
s# 计算 待定 参数 的 值 
D= MiD(1) =[]; 

D=D'， 

了 E=[- Ciones(1,n-1)]; 
C= inv(Ex 了 ) 关 马 # Di 
c=c3 
a=c(1);b=c(2)5 

s* 预测 后 续 数 据 
F=[];FG1) =RG); 

for i=2:(Cn+10) 


end 
G=[D;GG) =aRG); 
for i=2:Cn+10) 

GGi) =FCGi 一 PEGE-I1) 





207 234 


220.5256 270 
* 原 始 数据 累加 


多 生成 虹 加 矩阵 


FCGD = (RM1) - b/a)/exp(ax (i 一 1))+b/ay 


s# 得 到 预测 出 来 的 数据 


285]; 














end 

tl1 = 1995:2004; 
t2 = 1995:2014; 
Gia b 
plot(t1,a， o ,t2,G) 





当 输出 预测 值 、 发 展 系数 和 灰色 作用 量 
s# 原始 数据 与 预测 数据 的 比较 








550| 





200| 











ET 器 2006 2008 2010 2012 2014 


图 3-2 1995 一 2014 预计 长 江 排放 的 污水 量 


3.3.2 实例 二 :预测 与 会 代表 人 数 (CUMCM 2009 D) 


1， 问 题 描 述 


该 题 要 求 为 会 议 筹备 组 制订 一 个 预订 宾馆 客房 ,租借 会 议 室 .租用 客车 的 合理 方案 。 为 了 
解决 这 个 问题 ,需要 先 预 测 与 会 代表 的 人 数 。 预 测 的 依据 是 代表 回执 数量 及 往届 的 与 会 人 员 
数据 。 已 知 本 届 会 议 的 回执 情况 ( 见 表 3 - 电 及 以 往 几 届 会 议 代表 回执 和 与 会 情况 ( 见 表 3- 5)。 
要 解决 的 问题 是 :根据 这 些 数据 预测 本 届 与 会 代表 人 数 。 











表 3- 4 本 届 会 议 代表 的 回执 中 有 关 住 房 要 求 的 信息 单位 ;人 
要 求 合 住 1 合 住 2 合 住 3 独 住 1 独 住 2 独 住 3 
田 154 104 32 107 68 下 1 
78 48 17 59 28 19 





























说 明 : 表 头 第 一 行 中 的 数字 1.2、3 分 别 指 每 天 每 问 120 一 160 元 、161 一 200 元 .201 一 300 元 三 种 不 同 价格 的 房间 。 合 住 
是 指 要 求 两 人 合 住 一 间 。 独 住 是 指 可 安排 单 人 间或 一 人 单独 住 一 个 双人 间 。 


表 3-S5 以 往 几 届 会 议 代表 回执 和 与 会 情况 


单位 :人 





届 次 


第 一 届 


第 二 届 


第 四 届 





发 来 回执 的 代表 数 其 


315 


356 


408 


711 





发 来 回执 但 未 与 会 的 代表 数 基 


89 


115 


121 


213 





未 发 回执 而 与 会 的 代表 数 其 








57 





69 








104 











2. 问题 的 求解 


根据 表 3-4 的 数据 ,可 知 本 届 发 来 回执 的 数量 为 755。 根 据 表 3 - 5 的 数据 ,可 以 知道 发 


来 回执 但 未 与 会 的 代表 数 和 未 发 回执 但 与 会 的 代表 与 发 来 回执 数量 问 的 关系 。 


定义 1: 未 知 与 会 率 一 未 发 回执 但 与 会 的 代表 的 数量 /发 来 回执 的 代表 数量 


定义 2: 缺 席 率 = 发 来 回执 但 未 与 会 的 代表 数 /发 来 回执 的 代表 数量 
根据 以 上 定义 ,可 以 得 到 往届 的 缺席 率 和 未 知 与 会 率 , 如 表 3- 6 所 列 。 


表 3-6 往届 的 缺席 率 和 未 知 与 会 率 


























届 次 第 一 届 第 二 届 第 三 届 第 四 局 
缺席 率 0. 282540 0.323034 0. 296 569 0. 299 578 
未 知 与 会 率 0.180952 0. 193820 0.183824 0.146 273 





从 表 3- 6 可 以 看 出 ,缺席 率 一 直 保持 在 0. 3 左右 ,而 未 知 与 会 率 却 变 化 较 快 。 为 此 ,认为 
第 五 届 的 缺席 率 仍 为 0. 3。 这 样 ,缺席 的 人 数 为 755 人 X0.3 一 226.5 人 。 为 了 保守 起 见 , 对 


226.5 人 进行 向 下 取 整 , 即 缺席 的 人 数 为 226 人 。 


未 知 与 会 率 变化 相对 剧烈 ,不 适合 应 用 比例 方法 确定 ,同时 由 于 数据 有 限 , 所 以 应 用 灰色 
预测 方法 比较 合适 。 从 实际 问题 的 角度 ,认为 以 未 知 与 会 率 为 研究 对 象 较为 合适 。 将 往届 的 
未 知 与 会 率 数据 带 人 程序 P3 - 1, 并 对 输入 数据 和 预测 数据 做 相应 修改 ,可 很 快 得 到 本 届 的 未 
知 与 会 率 为 0. 1331 ,所 以 本 届 未 发 回执 但 与 会 的 代表 数量 为 755 人 X0.1331= 10054905 人 。 
同样 保守 考虑 ,向 上 取 整 为 101 人 。 这 样 就 可 以 预测 本 局 与 会 代表 的 数量 为 755 人 十 101 人 一 


226 大 =630 人 。 





第 4 章 
| 遗传 算法 及 其 MATLAB 实现 








4.1 遗传 算法 基本 原理 


4.1.1 人 工 智 能 算法 概述 


前 面 已 经 学 习 了 数学 建 模 中 常规 方法 的 使 用 .不 太 复 杂 问 题 的 MATLAB 规划 工具 箱 的 
求解 以 及 灰色 预测 理论 。 从 本 章 开始 ,将 陆续 接触 一 些 常 见 的 .相对 复杂 的 智能 算法 和 其 他 一 
些 高 级 方法 的 基本 理论 并 用 MATLAB 软件 实现 。 

海量 数据 、 信 息 高 速 公路 和 知识 爆炸 都 是 现代 社会 三 个 重要 的 基本 属性 ,数据 ,信息 和 知 
识 是 数字 时 代 的 关键 词 。 尽 管 知识 离 不 开 数据 和 信息 ,但 知识 与 数据 和 信息 并 不 是 同一 概念 。 
根据 美国 著名 管理 学 家 Richard L. Daft 的 经 典 著作 《组 织 理论 与 设计 ) 的 释义 ,数据 (data) 是 
简单 的 、 纯 粹 的 事实 和 数字 ,其 本 身 包含 的 使 用 价值 很 小 。 比 如 ,一 个 企业 可 能 有 项 数据 显示 
出 其 生产 的 某 种 产品 有 20 外 销 往 江苏 省 。 要 想 这 个 数据 有 用 ,就 必须 将 其 与 其 他 数据 横向 或 
纵向 联系 起 来 , 如 “这 批 产品 被 16 一 25 岁 的 青年 男女 买 走 ”, 将 其 处 理 为 信息 。 所 谓 信息 
Cinformation) 就 是 与 其 他 的 数据 联系 后 所 转化 成 的 具有 特定 用 途 的 具有 使 用 价值 的 数据 。 
知识 (knowledge) 就 更 进一步 了 , 它 是 指 将 某 三 信息 与 其 他 信息 相 联系 ,并 与 已 知 的 信息 进行 
比较 后 而 从 这 一 信息 中 得 出 的 结论 。 与 信息 和 数据 对 比 , 知 识 总 是 包含 着 人 或 者 智能 的 因素 。 
书本 能 承载 信息 ,但 是 这 些 必须 被 人 咀 收 并 付 诸 使 用 后 才能 成 为 知识 。 书 本 上 的 知识 是 建立 
在 已 有 的 信息 ` 亲 身体 验 、 直 党 和 理解 的 基础 上 的 , 它 包 含 着 所 掌握 的 信息 应 如 何 行动 的 认 知 
能 力 。 

智能 算法 是 智能 技术 领域 的 一 个 分 支 。 智 能 算法 出 现 的 原因 是 ,人 们 生活 在 知识 新 陈 代 
谢 速度 快 和 信息 繁杂 的 社会 里 ,需要 用 高 效 的 数据 挖掘 工具 从 自然 人 的 各 类 密码 ,账号 .消费 
习惯 所 形成 的 数据 中 提取 出 有 益 的 信息 和 知识 以 便 为 产品 的 销售 提供 决策 支持 ,以 及 从 企业 
的 各 种 海量 数据 报表 中 筛选 出 信息 和 知识 以 便 提 高 生产 效率 降低 成 本 ,或 者 用 计算 机 代替 人 
脑 进行 复杂 的 各 类 计算 。 如 何 从 海量 数据 中 采集 到 信息 ,又 如 何 从 信息 中 归纳 知识 ,是 一 项 高 
技术 并 且 劳 动 强度 很 大 的 工作 。 为 了 减轻 人 们 的 劳动 量 , 同 时 也 为 了 节约 社会 成 本 ,出现 了 一 
些 载 有 能 代替 人 类 脑力 劳动 的 算法 的 机 器 硬件 ,这 些 算法 被 称 为 智能 算法 ,它们 具有 一 些 类 似 
生命 体 智慧 的 特征 ,其 实 就 是 经 过 抽象 以 后 的 机 器 学 习 。 不 论 是 粒子 群 算法 还 是 蚁 群 算法 ,也 
不 论 是 遗传 算法 还 是 免疫 算法 ,更 不 论 是 神经 网 络 还 是 模拟 退火 算法 ,本 质 上 都 是 机 器 学 习 。 
机 器 学 习 具 有 一 个 显著 特征 一 一 机 械 性 。 吴 文俊 院士 在 (数学 机 械 化 》( 科 学 出 版 社 ,2003 年 3 
月 ) 一 书 中 提 到 ,数学 机 械 化 是 指数 学 在 运算 和 证 明 过 程 中 ,只 要 前 进 了 一 步 ,都 会 有 确定 的 、 
必然 的 下 一 步 ,直到 到 达 结论 ,整个 过 程 是 按照 既定 的 刻板 规律 进行 。 现 在 流行 的 用 计算 机 将 
解析 几何 转化 成 向 量 进行 证 明 就 是 20 世纪 70 一 80 年 代 吴 文俊 院士 首创 的 ,这 是 数学 机 械 化 
的 一 个 成 功 应 用 。 可 见 , 智 能 算法 、 机 器 学 习 和 数学 机 械 化 在 本 质 上 是 相通 的 ,可 以 被 看 成 同 














一 个 问题 的 不 同 侧面 。 

智能 算法 是 一 门 边缘 交叉 学 科 , 是 生物 .数学 等 多 学 科 的 完美 融合 。 各 学 科 交叉 之 处 往往 
展现 出 最 强劲 的 发 展 势头 。 思 格 斯 在 (自然 辩证 法 }》 里 就 说 过 :在 分 子 科学 和 原子 科学 的 接触 
点 上 ,物理 学 家 和 化 学 家 都 承认 自己 无 能 为 力 ,然而 应 当 在 这 点 上 期 待 最 大 的 成 果 。” 这 里 常常 
会 形成 一 门 新 的 科学 ,是 “新 科学 "的 “生长 点 ", 是 最 具 生命 力 ` 最 活跃 的 地 方 。 例 如 , 螨 烛 火 焰 
的 温度 总 在 略 带 淡 蓝 色 的 外 焰 温 度 最 高 , 焰 心 处 温度 最 低 ; 台 风 总 在 边缘 处 风力 最 大 , 风 眼 的 
风力 为 零 级 ;靠近 锅 边 的 水 总 是 最 先 沸腾 ;电荷 在 导体 的 尖端 边缘 的 地 方 容易 聚集 ;海岸 是 陆 
地 和 海洋 作用 最 强烈 的 地 方 ;磁铁 两 端的 NMS 极 总 是 磁性 最 强 。 

现代 的 智能 算法 已 经 很 多 了 ,不同 的 智能 算法 之 间 相 互 借鉴 ,不 断 融 合 ,差距 越 来 越 小 。 


4.1.2 站 传 算法 生物 学 基础 


在 一 定 的 时 间 里 ,有 一 群 免 子 , 其 中 一 些 比 另外 一 些 兔子 跑 得 快 ,而 且 更 聪明 ,这 些 兔子 被 
狐狸 吃 掉 的 可 能 性 比较 小 ,因此 它们 中 的 多 数 就 存活 下 来 并 繁殖 更 多 的 兔子 。 当 然 ,一 些 跑 得 
慢 而 愚蠢 的 兔子 也 会 存活 下 来 ,只 是 因为 它们 比较 侥幸 ,这 些 存活 的 兔子 群 开始 生育 。 生 育 的 
结果 是 兔子 遗传 材质 的 充分 融合 :一 些 跑 得 慢 的 兔子 生出 了 跑 得 快 的 兔子 ,一 些 跑 得 快 的 兔子 
生出 跑 得 更 快 的 ,一 些 聪明 的 兔子 生出 了 愚蠢 的 兔子 ,等 等 在 最 顶层 ,自然 界 不 时 地 变异 一 
些 兔 子 的 基因 材质 。 所 产生 的 小 兔子 平均 来 说 要 比 原始 的 群体 更 快 更 聪明 ,因为 从 狐狸 口中 
生存 下 来 的 父 代 多 数 是 跑 得 更 快 .更 聪明 的 兔子 。 同 样 ,狐狸 也 经 历 相似 的 过 程 ,否则 兔子 可 
能 跑 得 太 快 又 太 聪 明 以 致 狐狸 根本 抓 不 到 了 。 

兔子 的 生存 哲学 就 是 以 自然 选择 学 说 为 核心 的 现代 生物 进化 理论 ;其 基本 观点 是 :种 群 是 
生物 进化 的 基本 单位 ,生物 进化 的 实质 是 种 群 基因 频率 的 改变 。 基因 突 变 和 基因 重组 .自然 选 
择 及 隔离 是 物种 形成 过 程 的 三 个 基本 环节 ;通过 它们 的 综合 作用 ,种 群 产生 分 化 ,最 终 导致 新 
物种 的 形成 。 在 这 个 过 程 中 ,基因 突变 和 基因 重组 产生 生物 进化 的 原材料 ,自然 选择 使 种 群 的 
基因 频率 定向 改变 并 决定 生物 进化 的 方向 ,隔离 是 新 物种 形成 的 必要 条 件 。 

新 物种 形成 的 途径 和 方式 有 丙种 :渐变 式 和 爆发 式 。 渐 变 式 主要 通过 变异 的 逐渐 积累 而 
成 亚 种 ,再 由 亚 种 形成 一 侈 或 多 个 新 种 ,新 种 又 分 为 两 种 类 型 , 即 继承 式 新 种 形成 和 分 化 式 新 
种 形成 ;爆发 式 不 通过 亚 种 这 一 阶段 而 迅速 形成 新 的 物种 ,其 分 为 三 种 类 型 , 即 杂 交 产 生 新 种 ， 
染色 体 结构 变化 形成 新 种 和 多 倍 体 化 的 新 种 形式 。 遗 传 算法 杂交 了 渐变 式 和 爆发 式 的 两 种 
思想 。 


4.1.3 站 传 算法 的 实现 步骤 


1975 年 Michigan 大 学 的 J. Holland 教授 提出 “遗传 算法 "的 概念 并 出 版 了 颇 有 影响 的 著 
作 《Adaptation in Natural and Artificial Systems》。 遗 传 算法 是 模拟 达尔 文生 物 进 化 论 的 自 
然 选 择 和 备 德 尔 遗 传 学 机 理 的 生物 进化 过 程 的 计算 模型 ,是 一 种 通过 模拟 自然 进化 过 程 搜索 
最 优 解 的 方法 。J. Holland 教授 提出 的 遗传 算法 属于 简单 遗传 算法 (simple genetic algorithm 
,SGA) ,还 处 于 遗传 算法 发 展 进 程 的 初级 阶段 ,但 却 具 有 里 程 碑 式 的 意义 。 

SGA 由 编 解码 ,个体 适应 度 评估 和 遗传 运算 三 大 模块 构成 ,而 遗传 运算 又 包括 染色 体 复 
制 ,交叉 ,变异 甚至 倒 位 等 。 改 良 的 遗传 算法 和 融合 新 型 技术 的 遗传 算法 都 是 SGA 的 变异 形 
式 。 在 遗传 算法 中 ,定义 种 群 或 群体 为 所 有 编码 后 的 染色 体 集合 ,表征 每 个 个 体 的 是 其 相应 的 
染色 体 。 




















1. 编 码 
遗传 算法 的 编码 有 浮 点 编码 和 二 进 制 编码 两 种 ,这 里 只 介绍 二 进 制 编码 规则 。 二 进 制 编 
码 既 符合 计算 机 处 理 信息 的 原理 ,也 方便 了 对 染色 体 进 行 遗传 .编译 和 突变 等 操作 。 设 某 一 参 
数 的 取 值 范围 为 (L,U) ,使 用 长 度 为 上 的 二 进 制 编码 表示 该 参数 , 则 它 共 有 2 种 不 同 的 编码 。 
该 参数 编码 时 的 对 应 关系 为 
000000000000000000 王 0 一 三 
000000000000000001 王 1 一 工 十 9 
000000000000000010 王 2 一 工 十 26 
000000000000000011 王 3 一 世 十 36 
111111111111111111 一 2 一 1-U 
易 知 ; 
ia 
2. 解 码 
解码 的 目的 是 为 了 将 不 直观 的 二 进 制 数据 串 还 原 成 十 进 制 。 设 某 一 个 体 的 二 进 制 编码 为 


ZiOai-sst 访 , 则 对 应 的 解码 公式 为 
ET 寻 亚 王 
未 全 7 十 碎 25 
(252 = 


例如 ? 设 有 参数 xE[2, 和 4, 现 用 5 位 二 进 制 数 对 工 进行 编码 ,可 得 2 一 32 条 染色 体 ， 
00000,00001,00010，0001ly -00100，00101，00M0O，00TT1 
01000,01001,01010, :0IOUH ，01100，01101，01110，01111 
10000,10001,10010。10011，10100，10101，10110，10111 
11000,11001,110108 11011，111008 11101，11110，11111 
对 于 任意 三 进 制 数 据 串 只 要 代 大 译 码 公式 ,就 可 以 得 到 对 应 的 解码 ,如 到 :二 10101" 它 对 
应 的 十 进 制 值 为 





5 
了 02 二 1X2 十 IOX2 十 1X2 十 0X2 十 1X2% 一 21 


则 对 应 参数 的 值 为 
4 一 2 
条 一 1 
遗传 算法 的 编码 和 解码 在 宏观 上 可 以 对 应 生物 的 基因 型 和 表现 型 ,在 微观 上 可 以 对 应 
DNA 的 转录 和 翻译 两 个 过 程 。 
3. 交 配 01001011 10010101 
“交配 运算 是 使 用 单 点 或 多 点 进行 交叉 的 算 子 。 首先 用 随 
机 数 产生 一 个 或 多 个 交配 点 位 置 ,然后 两 个 个 体 在 交配 点 位 置 互 
换 部 分 基因 码 ,形成 两 个 子 个 体 。 例 如 ,有 两 条 染色 体 S; 王 0l00010f 1001D11 
01001011,S: 一 10010101 交换 其 后 4 位 基因 ,如 图 4- 1 所 示 。 
S 一 01000101,S' :一 10011011 可 以 被 看 做 是 原 染 色 体 S, 和 S。 
的 子 代 染 色 体 。 


工 一 2 十 21 X 





一 3.3548 


图 4-1 染色 体 基因 交配 范例 





恋 


“突变 运算 ”是 使 用 基本 位 进行 基因 突变 。 为 了 避免 在 算法 和 迭代 后 期 出 现 种 群 过 早 收敛 ， 
对 于 二 进 制 的 基因 码 组 成 的 个 体 种 群 ,实行 基因 码 的 小 几率 翻转 , 对 于 二 进 制 编码 即 0 变 为 
1, 而 1 变 为 0。 例 如 ,将 染色 体 S 王 11001101 第 3 位 上 的 0 变 为 1, 即 S 王 11001101 一 
11101101 一 S' 。S' 可 以 被 看 做 是 原 染色 体 S 的 子 代 染 色 体 。 

5. 倒 位 

除了 交配 和 突变 之 外 ,对 于 复杂 的 问题 可 能 需要 用 到 * 倒 位 ”, 其 对 应 的 运算 亦 被 称 为 “ 倒 
位 运算 ”。 倒 位 是 指 一 个 染色 体 某 区 段 正常 排列 顺序 发 生 180" 的 颠倒 ,造成 染色 体内 的 DNA 
序列 重新 排列 , 它 包 括 臂 内 倒 位 和 臂 间 倒 位 。 倒 位 纯 合体 不 影响 个 体 的 生活 力 , 只 是 改变 了 染 
色 体 上 的 相 邻 基因 位 置 , 从 而 某 些 表现 型 发 生 位 置 效应 ,同时 也 改变 了 与 相 邻 基因 的 交换 值 。 
倒 位 杂 合 体 则 不 然 ,其 生育 力 降低 。 染 色 体 上 的 区 段 可 能 一 次 又 一 次 发 生 倒 位 ,上 且 通 过 自 交 出 
现 不 同 的 倒 位 纯 合体 ,致使 它们 与 其 原来 的 物种 不 能 交配 ,形成 生殖 隔离 ,结果 产生 新 族群 或 
变种 。 

倒 位 运算 是 与 倒 位 概念 相似 的 运算 规则 。 例 如 ,染色 体 S 就 是 S 经 过 倒 位 运算 以 后 得 到 
的 染色 体 基 因 编 码 , 其 中 倒 位 的 部 分 即 S 中 标记 下 画 线 的 部 分 。 








S = 1001011011101110011010101001=S' 一 100101Hi00lo11001110111101001 
6. 个 体 适应 度 评估 
自然 界 中 能 够 适应 环境 的 生物 有 更 多 的 机 会 存活 下 来 ;这 种 筛选 机 制 类 似 图 4 二 2 所 示 现 
象 :在 由 正六 边 形 搭建 的 三 角形 区 域 的 顶部 投掷 一 T 


些 光滑 的 木 块 ,这 些 木 块 经 由 白色 的 缝 辽 坠落 底部 ， 
显然 落 在 底部 中 间 的 木 块 要 比 落 在 两 端的 木 块 多 ， 
因为 木 块 有 更 多 的 路 径 坠 落 在 底部 的 中 间 区 域 , 所 
以 有 更 大 的 几率 落 在 中 间 。 落 在 各 个 区 域 的 几率 对 
应 遗传 算法 中 各 条 染色 体 被 遗传 到 下 一 代 的 几率 ， 
其 瞧 落 的 位 置 对 应 自 变 量 取 值 。 遗 传 算法 依照 与 个 
体 适应 度 成 正比 的 几率 决定 当前 种 群 中 各 个 个 体 遗 
传 到 下 一 代 群 体 中 的 机 会 个 体 适 应 度 大 的 个 体 更 
容易 被 遗传 到 下 一 代 。 通 常 , 求 目标 函数 最 大 值 的 
问题 可 以 直接 把 目标 函数 作为 检测 个 体 适 应 度 大 小 
的 函数 。 

7. 复 制 

复制 运算 是 根据 个 体 适 应 度 大 小 决定 其 下 代 遗 传 的 可 能 性 。 若 设 种 群 中 个 体 总 数 为 N， 
个 体 ; 的 适应 度 为 矿 , 则 个 体 i 被 选取 的 几率 
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图 4-2 正六 边 形 筛选 机 制 


1 
和 

当 个 体 复制 的 几率 决定 后 ,再 产生 [0， 可 区 辣 的 均匀 随机 数 来 决定 哪些 个 体 参加 交 醒 。 若 
个 体 适 应 度 高 , 则 被 选取 的 几率 已 就 大 , 则 可 能 被 多 次 选中 , 它 的 遗传 基因 就 会 在 种 群 中 扩 
散 ;: 若 个 体 的 复制 几率 小 , 则 会 被 逐渐 淘汰 .。 

下 面 以 一 个 经 典 案例 具体 说 明 遗 传 算法 的 作业 过 程 。 


书 








【 例 4-1] 求 maxFCzi,zrz) 一 21.5 二 risin(4xri) 十 zasin(20rzrs) 
| 一 3.0 私 屏 委 12.1 
141 委 委 5.8 

函数 fCri ,zs ) 的 三 维 图 形 如 图 4- 3 所 示 。 


光志 





图 4-3 函数 Jrivza) 的 三 维 图 形 


(1) 编码 

首先 要 进行 编码 工作 ,即将 变量 转换 成 二 进 制 数 串 。 数 串 的 长 度 取决 于 所 要 求 的 精度 。 
例如 ,变量 z 的 区 间 是 (L,UD) ,要 求 的 精度 是 小 数 点 后 4 位 ,也 就 意味 着 每 个 变量 应 该 被 分 成 
至 少 (LAL 门 X10: 企 部 分 。 对 一 企 变 量 的 二 进 制 数 串 位 数 用 以 下 公式 计算 ; 

2 SU 荔 罗 OF 2 一 1 

本 例 精度 要 求 保留 小 数 点 后 4 位 , 则 目标 函数 的 两 个 自 变量 z: 及 zs 所 构成 的 染色 体 数 
串 可 以 表示 如 下 ， 
一 3.0 委 委 12.1 
[12.1 一 (一 3.0)]X10: 一 151000 
2 <151000 委 2 一 1 
ml 一 18 
4.1 和 mm 入 5.8 
(5.8 一 4.1) X 10: 一 17000 
2 天 17000 委 2 一 1 
ma 一 15 
六 一 ml 十 ma 一 18 十 15 一 33 
本 例 中 任 一 染色 体 数 串 都 是 33 位 , 即 000001010100101001101111011111110。 
以 上 编码 前 18 位 表示 zi, 后 15 位 表示 立 , 如 表 4-1 所 列 。 





表 4-1 染色 体 编码 











自 变 基 二 进 制 十 进 制 实际 值 
羡 000001010100101001 5417 一 2.687969 
三 101111011111110 24318 5.361653 




















则 二 进 制 转化 成 十 进 制 为 
了 


一 3.0+5417X 世 1 
辣 一 和 1 十 24318X 5 人 1 
假设 初始 种 群 中 有 10 个 个 体 ,其 染色 体 可 随机 生成 如 下 : 

U = [000001010100101001 101111011111110] 
Us = [ool110101110011000 000010101001000] 
Us = [111000111000001000 010101001000110] 
U, = [1o0110110100101101 000000010113001] 
Us = [000010111101100010 001110001401000] 
Us = [111110101011011000 0000101100ilo01] 
U; = Clolooololl11l0001 oo1100111011101] 
柄 二 [oolol1010100001100 010Tlo011001100] 
夸 =[ullloool011101100 0Tilol000111101] 
Un = [lll1010011101010Io o00010101101010J 
相对 应 的 十 进 制 的 实际 值 Lzi ,zs 为 
U = [一 2. 687969,5.361653]， Us = [0.474101,4. 17014 和 ] 
Us = [10. 419457,4:661861]， 。 U, 一 虑 6. 159951,4. 109 598] 
项 = 章 o 了 js6.772o 柯 用 JEIDi7Ros4 ,人 5 
疝 芋 [0842067:5i21702]，， Us 户 [0.330256,4.694979 
Us = [11. 677267,4.873501]， 。 Uns 一 [11.446 273,4. 171908] 
(2) 评价 个 体 造 应 度 
对 一 个 染色 体 数 串 的 适应 度 的 评价 由 下 列 三 个 步 又 组 成 。 
@ 将 染色 体 串 进行 反 编 码 (解码 ) ,转换 成 真实 值 , 即 
愉 一) 一 123 
Q@ 评价 目标 函数 Fr ) 。 
@@ 将 目标 函数 值 转 为 适应 度 。 对 于 极 大 值 问题 ,适应 度 可 作为 目标 函数 值 : 
evaldD = FaD， 大 一 123 
在 遗传 算法 中 ,评价 函数 扮演 自然 进化 中 环境 的 角色 , 它 通 过 染色 体 的 适应 度 对 其 进行 评 
价 。 上 述 染 色 体 的 适应 度 值 如 下 : 
eval(U) 一 FC- 2. 687969,5. 361653) 一 19. 805119 
evalUa) 一 (0.474101,4.170144) 一 17. 370896 
eval(Ua) 一 /10. 419457,4. 661461) 一 9. 590546 
eval(U) 一 76. 159951,4.109598) 一 29.406122 


.0) 





一 一 2.687969 


一 5.361653 





eval(Us) 一 太一 2. 301286,4.477 282) 一 15. 686 091 
eval(Us) 一 j(11.788084,4.174346) 一 11. 900 541 
eval(U7?) 一 8 (9. 342067,5.121702) 一 17. 958717 
eval(Us) 一 太一 0.330256,4. 694977) 一 19. 763 190 
eval(U。) = j(11. 671267,4. 873 501) 一 26. 401 669 
eval(Un) 一 (11. 446 273,4. 171 908) 一 10. 252480 

依照 染色 体 的 适应 度 值 进行 新 种 群 的 复制 ,步骤 如 下 : 

@ 计算 染色 体 LT 的 适应 度 值 

eval(Us) 一 fx， 大 一 1 2 
@ 计算 种 群 的 适应 度 值 总 和 


ps 


忆 一 Zeval(U) 
@@ 计算 每 个 染色 体 被 复制 的 概率 


己 二 evaldw 
后 
@@ 计算 每 个 染色 体 被 复制 的 累积 概率 
QE 已 P 
所 


(3) 新 种 群 复制 
依照 轮 盘 选择 法 ,转动 轮 盘 10 次 (种 群 中 有 .10 条 染色 体 ), 每 次 选择 一 个 作为 新 种 群 的 染 
色 体 。 假 设 10 次 中 产生 的 0 一 1 随机 数 序列 如 下 : 
0.301431 0.322062 :0.766503 0.881893 0.350871 
0.538392 0.177618 0.343242 0.032685 0.197577 
根据 以 上 的 计算 方法 ,可 以 先 计算 出 种 群 中 每 个 染色 体 的 适应 度 和 概率 ,如 表 4- 2 所 列 。 
利用 计算 机 模拟 轮 盘 选择 法 ,假设 计算 机 产生 10 个 [0, 蕊 区间 的 随机 数列 如 下 : 
表 4 -2 .种 群 每 条 染色 体 的 适应 度 、 被 复制 概率 和 被 复制 的 累积 概率 



































染色 体 适应 度 值 书 Q， 
Un 19.805119 0.111180 0.111180 
Uz 17. 370896 0.097515 0.208695 
Ua 9.590 546 0.053839 0. 262534 
Ua 29. 406 122 0.165077 0.427611 
Us 15. 686091 0.088057 0.515668 
Us 11. 900 541 0.066 806 0. 582475 
U7y 17.958717 0.100815 0.683 290 
Us 19.763 190 0.110945 0.794234 
Us 26. 401669 0.148211 0.942446 
Une 10. 252 480 0.057 554 1.000000 




















第 1 个 随机 数 为 0. 301 431, 大 于 Q: 小 于 Q, ,所 以 Us 被 选中 ; 
第 2 个 随机 数 为 0. 322 062, 大 于 Q, 小 于 Q, ,所 以 U 再 次 被 选中 ; 


第 3 个 随机 数 为 0. 766 503, 大 于 Q; 小 于 Qs ,所 以 Us 被 选中 ; 
第 10 个 随机 数 为 0. 197 577, 大 于 Qi 小 于 Q: ,所 以 Us 被 选中 ; 
依照 轮 盘 选择 法 ,新 种 群 的 染色 体 组 成 如 下 : 
Un = [100110110100101101 000000010111001](U4) 
Us = [100110110100101101 000000010111001](U4) 
Us = [001011010100001100 010110011001100](Us) 
Us = [111110001011101100 011101000111101]GU,) 
Us = [100110110100101101 000000010111001]CU4) 
Us = [110100010111110001 001100111011101]CU;，) 
U; = [001110101110011000 000010101001000](U:) 
Us = [100110110100101101 000000010111001](U:) 
Un = [000001010100101001 101111011111110](U,) 
Un = [001110101110011000 000010101001000](U2: ) 
这 种 轮 盘 选择 法 的 机 理 是 :染色 体 的 适应 度 大 意味 着 [QQ5i] 区 间 路 度 就 大 ,随机 数 发 
生 器 产生 的 均匀 随机 数 就 会 有 更 大 的 概率 落 在 较 大 长 度 的 [QsQ,+,] 区 间 里 ,这 样 具 有 较 大 
已 值 的 染色 体 自 然 更 有 机 会 复制 到 下 一 代 。 
(4) 新 种 群 交配 
O 交配 染色 体 数量 的 确定 
交配 染色 体 的 数量 等 于 染色 体 总 量 乘 以 交配 概率 。 这 里 假设 交配 概率 己 . 为 0. 25 ,染色 
体 总 量 为 10 条 ,所 以 参加 交配 的 染色 体 数 量 为 [2. 5] 条 。 符 号 [ ] 表 示 取 整 , 这 里 取 整 数 2, 即 
交配 的 染色 体 数目 为 2 条 。 
四 交配 染色 体 对 象 的 确立 
用 计算 机 产生 [0, 匡 区 间 的 15 个 随机 数 如 下 : 
0.625724 只 0.266823， 0.288644 丰 0.295114。， 0.163274 
0.567461 “707085940 0.392865 ”0.770714 。 0.548656 
假定 其 分 别 对 应 由 一 Duo 这 10 个 个 体 , 则 其 中 低 于 交配 概率 0. 25 的 Us 和 LU， 参加 交配 。 
这 样 操作 的 原因 是 :交配 概率 越 低 , 低 于 交配 概率 以 下 的 随机 数 的 数量 就 越 少 , 所 以 参加 交配 
的 染色 体 数量 与 交配 概率 可 能 会 成 正比 。 
轩 在 交配 池 发 生 交 配 
染色 体 Us 和 LU 被 选中 作为 交配 的 父辈 ,交配 点 的 选择 以 随机 数 产生 。 交 配 的 种 类 有 单 
点 交配 和 多 点 交配 ,这 里 取 单 点 交配 。 计 算 机 随机 生成 一 个 介 于 0 一 32 的 整数 (因为 整个 染色 
体 数 串 的 长 度 为 33) 。 假 设 所 产生 的 整数 为 1 ,那么 两 个 染色 体 自 1 位 置 ( 即 二 进 制 串 的 第 二 
位 ) 开 始 分 割 ,在 染色 体 1 位 置 右 端 部 分 进行 交换 而 生成 新 的 子 辈 染色 体 , 即 
Us = [1 0011 0110 1001 0110 1000 0000 1011 1001] 
U; = [0 0111 0101 1100 1100 0000 0101 0100 1000] 
U3 = [1 0111 0101 1100 1100 0000 0101 0100 1000] 
U7 = [0 0011 0110 1001 0110 1000 0000 1011 1001] 




















(5) 基因 突变 


依照 突变 运算 规则 并 假设 突变 几率 P。 为 0. 01, 亦 即 种 群 内 所 有 基因 都 有 0. 01 的 概率 进 
行 突变 。 在 本 例 中 共有 33X10 一 330 个 基因 , 即 希望 每 一 代 中 有 3. 3 个 突变 基因 ,每 个 基因 的 


突变 几率 是 均等 的 。 因 此 ,将 产生 330 个 介 于 0 一 1 之 间 的 随机 数 ( 需 编 号 ) ,然后 


将 该 随机 数 


小 于 0.01 者 选 出 ,并 将 其 对 应 的 基因 值 加 以 翻转 ,假设 330 次 中 产生 0 一 1 之 间 的 随机 数 ,其 


值 小 于 0.01 者 如 表 4- 3 所 列 。 


表 4-3 基因 突变 位 置 



































染色 体位 村 基因 位 数 随机 数 
表 6 0..009857 
5 32 0.003113 
7…1 生 0.000946 
329 十 33 一 10…32 10 32 0.001282 








表 4-3 第 1 列 显示 的 是 具体 在 哪些 染色 体 以 及 在 染色 体 的 什么 位 置 进行 了 突变 。 例 如 
第 1 行 "105 二 33 一 4…6” 表 示 的 是 在 第 4 条 染色 体 第 6 个 基因 生发 生 了 突变 ,因为 第 4 条 染色 
体 第 6 个 基因 对 应 的 基因 编号 是 33X (4 一 1) 十 6 一 105。 


在 突变 后 ,最 终 新 种 群 的 染色 体 组 成 如 下 : 


琴 三 [loo110110100101101000900010111001] 
灰 三 [lo0l110110100101101 000000010111001] 
Us = [001011010100001100 010110011001100] 
U, = [1111110010I101100 011101000111101] 
Us = [loll1010il10011000 000010101001010] 
Us = [lloioo0l0111110001 00 划 00111011101] 
项 三 LIoOiiOiioW0ioll01 0o000001011100 切 
US Do0llo1101o01011o1 00000001011100 刁 
Us 关 -Lo00001010100101001 101111011111110] 
Un = [001110101110011000 000010101001010] 


新 一 代 的 相对 应 实际 值 Lri,z*] 和 适应 度 值 如 下 : 


eval(Ui ) 一 6. 159 951,4. 109 598) 一 29. 406 122 
eval(Ua ) 一 8(6. 159 951,4. 109 598) 一 29. 406 122 


eval(Us) 一 F( 一 0.330256,4. 694977) 一 19.763 190 


eval(Uas ) 一 j(11. 907 206,4. 873 501) = 5. 702781 
eval(Us) 一 (8.024130,4. 170248) 一 19. 910 251 
eval(Ue) 一 帮 (9. 342067,5. 117020) 一 17. 958717 
eval(U; ) 一 8(6. 159 951,4. 109 598) 一 29. 406 122 
eval(Us ) 一 8(6. 159 951,4. 109 598) 一 29. 406 122 


eval(Us) 一 帮 一 2. 687969,5. 361653) 一 19. 805 119 
eval(Un) 一 帮 0.474101,4. 170248) 一 17. 370 896 
至 此 ,已 完成 遗传 算法 的 第 一 代 流 程 。 依 此 选 代 , 在 第 451 代 得 到 对 应 最 大 目标 函数 值 的 





Us = [111110000000111000 111101001010110] 

相应 实际 值 Lriyz?] 一 [11. 631 407,5. 724 824] ,适应 度 值 eval(U) 一 沁 (11. 631 407,5. 724 824) 一 
38.818 208 

从 以 上 实例 可 以 得 出 两 点 结论 : 

四 遗传 算法 本 质 上 是 一 种 启发 式 的 随机 搜索 算法 ,所 以 由 遗传 算法 得 出 的 结果 每 次 都 不 
尽 相 同 。 

@ 自 变量 在 给 定 的 约束 条 件 下 进行 了 无 缝 编码 ( 即 这 种 编码 方式 能 够 表达 解 空间 中 的 所 
有 可 行 解 ), 所 以 从 理论 上 讲 , 遗 传 算法 总 有 很 多 机 会 得 到 全 局 最 优 结果 而 不 是 局 部 最 优 结果 。 


4.1.4 遗 传 算法 的 拓展 


1. 协同 进化 遗传 算法 

遗传 算法 是 建立 在 生物 进化 理论 和 染色 体 遗 传 变异 基础 之 上 的 ,该 理论 遵循 优胜 劣 汰 的 
自然 选择 机 制 ,本 质 是 “ 胜 者 为 王 , 败 者 为 宅 ? 的 生物 种 群 内 部 激烈 的 竞争 关系 。 但 实际 上 , 生 
物 在 进化 过 程 中 ,除了 竞争 关系 ,还 有 协作 ,寄生 等 关系 。 生 物 与 生物 之 间 生物 与 环境 之 间 存 
在 着 “相似 相 溶 和 相反 相 成 ”的 和 谐 共 处 关系 。 所 以 * 基 于 达尔 文 进化 论 和 备 德 尔 遗 传 定律 的 
遗传 算法 不 可 和 避免 地 存在 着 弱点 ,这 为 标准 遗传 算法 处 理 异 常 复杂 的 优化 问题 时 显得 无 能 为 
力 得 到 优 证 。 

协同 进化 的 概念 最 早 由 Ehrlich 和 Raven 讨论 植物 和 植 食 昆虫 (蝴蝶 7 相互 之 间 的 进化 影 
响 时 提出 的 。 借 鉴 生态 系统 中 的 协同 和 协作 关系 ,Jason 提出 了 协同 进化 计算 机 制 , 并 给 出 了 
系统 进化 算法 的 具体 数学 模型 。 相 关 协 同 进化 理论 和 MATLAB 源 程序 可 参阅 文献 [4 。 

2. 遗传 算法 与 弗 洛 伊 德 梦 的 解析 法 

一 般 认 为 梦 是 这 样 形 成 的 :白天 ,大 的 左 大 脑 半球 工作 最 忙 , 它 从 感觉 器 官 中 接受 信息 , 进 
行 加 工 ,解决 可 能 出 现 的 各 种 问题 之 同时 , 右 大 脑 半 球 也 悄悄 地 工作 ,神经 元 记录 下 那些 自己 
来 不 及 考虑 的 情绪 和 信息 .。 当 大 们 熟睡 后 ,切换 到 右 大 脑 “ 主 持 ? 工 作 * 白 天 印 在 有 半球 里 的 全 
部 潜意识 以 虚拟 梦 的 形式 进行 释放 。 

如 果 将 梦 的 成 形 过 程 进行 离散 , 梦 由 以 下 程序 组 建 而 成 : 梦 的 改装 、 梦 的 置换 . 梦 的 替代 、 
梦 的 凝 缩 . 梦 的 转移 、 梦 的 表现 、 梦 的 仿 同 、. 梦 的 颠倒 和 梦 的 选择 ,经 由 这 些 步 又 对 梦 思 进 行 编 
译 ,编译 后 的 梦 内 容 一 般 难以 识别 ,而 后 通过 大 脑 的 整合 作用 将 这 些 代码 碎片 进行 链接 ,形成 
连贯 的 梦 的 情节 ,这 就 是 基于 弗 洛 伊 德 思维 的 梦 的 建 模 过 程 。 译 梦 则 是 梦 的 成 形 过 程 的 道 
运算 。 

着 重 简介 一 下 “ 梦 的 仿 同 ”成形 原 理 。 梦 的 仿 同 是 指 将 梦 思 中 的 具有 某 一 相同 属性 的 人 或 
物 的 性 状 都 集 于 某 一 人 或 某 一 物 上 ,这 个 过 程 称 为 仿 同 (集锦 )。 仿 同 是 用 在 人 身上 ,而 集锦 则 
是 对 事物 的 统一 。 不 过 集锦 亦 可 施用 在 人 身上 。 在 仿 同 作用 里 ,只 有 和 共同 元 素 相连 的 人 才 
能 表现 在 梦 的 显 意 里 ,其 他 人 则 被 抑制 了 。 这 跟 孟 德尔 的 遗传 定律 原理 是 一 样 的 , 当 显 性 基因 
存在 的 时 候 , 隐 性 基因 就 不 能 表现 出 隐 性 性 状 ,除非 在 完全 隐 性 的 环境 里 。 梦 在 形成 过 程 中 都 
要 经 过 审查 制度 。 一 些 有 违 共识 和 有 损 伦 理 常 岗 的 人 和 物 的 属性 是 逃 不 过 梦 的 审查 的 。 尽 管 
人 在 潜意识 里 都 有 兽性 的 一 面 ,在 利 比 多 的 驱使 下 都 有 原始 冲动 的 倾向 ,但 是 人 类 在 后 天 的 学 
习 过 程 中 受到 了 文化 舆论、 礼仪. 法律 .传统 等 方面 的 教化 和 约束 ,在 实际 生活 中 形成 了 牢固 
的 自我 约束 意识 ,即便 在 梦 中 这 种 意识 也 没有 消除 ,表现 为 梦 的 审查 制度 。 在 梦 中 如 果 不 能 逃 











脱 梦 的 审查 制度 的 素材 演变 成 隐 性 ,与 另 一 个 能 够 逃脱 审查 制度 的 人 或 物 集成 在 一 起 ,如 此 便 
能 逃脱 审查 ,这 就 是 仿 同 (集锦 ) 在 梦 形成 中 所 起 的 作用 。 梦 在 形成 过 程 中 ,就 这 样 伴随 着 审查 
与 欺骗 的 相互 博弈 。 

弗 洛 伊 德 梦 的 解析 法 虽然 隶属 于 心理 学 ,遗传 算法 是 人 工 智能 技术 的 一 个 分 支 ,表面 看 起 
来 两 者 之 间 没 有 任何 耦合 关联 ,但 是 仔细 探究 可 知 二 者 有 众多 相似 之 处 。 例 如 , 梦 的 置换 和 替 
代 可 以 类 比 于 遗传 算法 中 染色 体 上 基因 片段 的 交叉 和 变异 ; 梦 的 表现 可 以 类 比 遗 传 算法 中 生 
物 表现 型 向 基因 型 映射 的 过 程 , 也 就 是 遗传 算法 中 种 群 编码 的 过 程 ; 梦 的 颠倒 可 以 类 比 遗 传 算 
法 中 倒 位 运算 …… 其 实 梦 不 仅 与 遗传 算法 关系 密切 ,如 果 是 基于 弗 洛 伊 德 的 解梦 理论 , 梦 还 与 
小 波 分 析 、 神 经 网 络 . 人 工 免 疫 算法 等 均 有 联系 。 在 此 将 弗 洛 伊 德 梦 的 解析 法 与 遗传 算法 联系 
起 来 ,并 非 牵强 附会 ,而 是 希望 能 起 到 抛砖引玉 .举一反三 的 作用 ,希望 能 在 梦 与 人 工 智能 技术 
之 间 找 到 契合 点 ,开发 出 一 种 “算法 仿 同 ” 的 高 性 能 算法 ,希望 这 种 算法 有 朝 一 日 能 “愿望 达 
成 ”。 事 实 上 ,不 同 智能 算法 之 间 已 经 相互 仿 同 .相互 借鉴 相互 融合 ,有 异曲同工 的 趋势 ,差距 
越 来 越 小 了 。 











4.2 遗传 算法 的 MATLAB 程序 设计 


4.2.1 程序 设计 流程 及 参数 选取 


1， 遗传 算法 的 程序 设计 伪 代 码 

BEGIN 
t= 0 和 遗传 代数 
初始 化 PC(t) 初 始 化 种 群 或 染色 体 
计算 P(t) 的 适应 值 ; 


while (不 满足 停止 准则 ) do 
begin 
于 下 
估 PCt 本 ) 电 选择 PCDNi fs 选择 
重组 P(t) * 交叉 和 变异 
计算 P(t) 的 适应 值 
end 
了 END 


2 遗传 算法 的 参数 设计 原则 

在 单纯 的 遗传 算法 中 ,也 并 不 总 是 收敛, 即使 在 单 峰 或 单调 也 是 如 此 。 这 是 因为 种 群 的 进 
化 能 力 已 经 基本 丧失 ,种群 早 熟 。 为 了 避免 种 群 的 早熟 ,参数 的 设计 一 般 遵从 以 下 原则 

(1) 种 群 的 规模 

当 群 体 规模 太 小 时 ,很 明显 会 出 现 近亲 交配 ,产生 病态 基因 ,而 且 造成 有 效 等 位 基因 先天 
缺乏 ,即使 采用 较 大 概率 的 变异 算 子 ,生成 具有 竞争 力 高 阶 模式 的 可 能 性 仍 很 小 ,况且 大 概率 
变异 算 子 对 已 有 模式 的 破坏 作用 极 大 。 同 时 遗传 算 子 存在 随机 误差 (模式 采样 误差 ) ,妨碍 小 
群体 中 有 效 模式 的 正确 传播 ,使 得 种 群 进 化 不 能 按照 模式 定理 产生 所 预测 的 期 望 数量 :种群 规 
模 太 大 ,结果 难以 收敛 且 浪费 资源 ,稳健 性 下 降 。 种 群 规模 的 一 个 建议 值 为 0 一 100。 


(2) 变异 概率 

当 变 异 概率 太 小 时 ,种 群 的 多 样 性 下 降 太 快 ,容易 导致 有 效 基因 的 迅速 丢失 且 不 容易 修 
补 ; 当 变异 概率 太 大 时 ,尽管 种 群 的 多 样 性 可 以 得 到 保证 ,但 是 高 阶 模式 被 破坏 的 概率 也 随 之 
增 大 。 变 异 概率 一 般 取 0. 0001 一 0. 2。 

(3) 交配 概率 

交配 是 生成 新 种 群 最 重要 的 手段 。 与 变异 概率 类 似 ,交配 概率 太 大 容易 破坏 已 有 的 有 利 
模式 ,随机 性 增 大 ,容易 错失 最 优 个 体 ;交配 概率 太 小 不 能 有 效 更 新 种 群 。 交 配 概率 一 般 取 
0.4 一 0.99。 

(4) 进化 代数 

进化 代数 太 小 ,算法 不 容易 收敛 ,种 群 还 没有 成 熟 :代数 太 大 ,算法 已 经 熟练 或 者 种 群 过 于 早熟 
不 可 能 再 收敛 ,继续 进化 没有 意义 ,只 会 增加 时 间 开支 和 资源 浪费 。 进 化 代数 一 般 取 100 一 500。 

(5) 种 群 初 始 化 

初始 种 群 的 生成 是 随机 的 。 在 初始 种 群 的 赋予 之 前 ,尽量 进行 一 个 大 概 的 区 间 估计 ,以 免 
初始 种 群 分 布 在 远离 全 局 最 优 解 的 编码 空间 ,导致 遗传 算法 的 搜索 范围 受到 限制 ,同时 也 为 算 
法 减轻 负担 。 

3. 适应 度 函数 的 调整 

(1) 在 遗传 算法 运行 的 初期 阶段 

群体 中 可 能 会 有 少数 几 个 个 体 的 适应 度 相 对 其 他 个 体 来 说 非常 高 。 落 按照 常用 的 比例 选 
择 算 子 确定 个 体 的 遗传 数量 时 , 则 这 几 个 相对 较 好 的 个 体 将 在 下 一 代 群 体 让 可 有 很 高 的 比例 ， 
在 极端 情况 下 或 当 群 体 规模 较 小 时 ,新 的 群体 甚至 完全 由 这 样 的 少数 几 个 个 体 组 成 。 这 时 交 
配 运算 就 起 不 了 什么 作用 ,因为 相同 的 两 个 个 体 不 论 在 何 处 发 生 交叉 行为 都 永远 不 会 产生 新 
的 个 体 。 这 样 就 会 使 群体 的 多 样 性 降低 ,容易 导致 遗传 算法 发 生 早熟 现象 (或 称 早期 收敛) ,使 
址 传 算法 所 求 的 解 停留 在 某 一 局 部 最 优点 正六 因此 ;希望 在 遗传 算法 运行 的 初期 阶段 ,算法 能 
够 对 一 些 适 应 度 较 高 的 个 体 进行 控制 ,降低 其 适应 度 与 其 他 个 体 适 应 度 之 问 的 差异 程度 ,从 而 
限制 其 复制 数量 ,以 维护 群体 的 多 样 性 。 

(2) 在 遗传 算法 运行 的 后 期 阶段 

群体 中 所 有 个 体 的 平均 适应 度 可 能 会 接近 于 群体 中 最 佳 个 体 的 适应 度 。 也 就 是 说 ,大 部 
分 个 体 的 适应 度 和 最 佳 个 体 的 适应 度 差 异 不 大 ,它们 之 间 无 竞争 力 ,都 会 有 以 相 接近 的 概率 被 
遗传 到 下 一 代 的 可 能 性 ,从 而 使 得 进化 过 程 无 竞争 性 可 言 , 只 是 一 种 随机 的 选择 过 程 。 这 将 导 
致 无 法 对 某 些 重点 区 域 进 行 重点 搜索 ,从 而 影响 遗传 算法 的 运行 效率 。 因 此 ,希望 在 遗传 算法 
运行 的 后 期 阶段 ,算法 能 够 对 个 体 的 适应 度 进行 适当 的 放大 ,扩大 最 佳 个 体 适 应 度 与 其 他 个 体 
适应 度 之 间 的 差异 程度 ,以 提高 个 体 之 间 的 竞争 性 。 


4.2.2 MATLAB 站 传 算法 工具 箱 


1，GA 工具 箱 版 本 
为 了 省 略 艰深 难 懂 的 遗传 算法 数学 理论 和 降低 程序 开发 的 难度 ,MATLAB 软件 已 经 将 
遗传 算法 命令 进行 了 封装 ,做 成 专门 的 遗传 算法 工具 箱 一 GA Toolbox, 方 便 用 户 调用 。 关 
于 遗传 算法 工具 箱 ,需要 说 明 三 点 : 





